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ABSTRACT 

There are strong evidence for powerful jets in the low/hard state of black-hole X-ray binaries 
(BHXRBs). Here, we present a model in which electrons are accelerated once at the base of the 
jet, and are cooled by synchrotron emission and possible adiabatic energy losses. The accelerated 
electrons assume a Maxwellian distribution at low energies and possible energetic power law tail. 
These assumptions yield to a wealth of spectra, which we study in details. We identify critical values 
of the magnetic field, and five transition frequencies in the spectra. In particular, we show that: (I) 
the decay of the magnetic field along the jet enables, for wide jets, production of flat radio spectra 
without the need for electrons re-acceleration along the jet. (II) An increase of the magnetic field 
above a critical value of ~ 10^ G leads to a sharp decrease in the flux at the radio band, while the flux 
at higher frequencies saturates to a constant value. (Ill) For strong magnetic field, the flux decays 
in the optical/UV band as F^, oc v~^/'^^ irrespective of the electrons initial distribution. (IV) For 
-Bq ~ 10"* G, the X-ray flux gradually steepens. (V) With adiabatic energy losses, flat spectrum can 
be obtained only at a limited frequency range, and under certain conditions (VI) For narrow jets, 
r{x) oc with a < f/2, flat radio spectrum cannot be obtained. We provide full description of 
the spectrum in the different scenarios, and show that our model is consistent with the key observed 
properties of BHXRBs. 

Subject headings: plasmas — radiation mcchanisms:non-thcrmal — stars:winds, outflows — X- 
rays: binaries 



1. INTRODUCTION 

Observations of the low-hard state of many black hole 
X-ray binaries (BHXRBs) reveal a rich broadband spec- 
trum, extending from the radio to the hard X-rays (for 
reviews, see, e . g..lLiang Ill998t iMcClintock fc Remillard I 
120061 : IVrtilek I |2008D. Typically, the radio spectrum in 
the GHz - THz range is fla t (F^ oc i/"" with a ^ 0) 
or sli ghtly inverted {a < 0) ijHvnes et al. 1 120001 : iFender I 
I2OOII [2006) . This s pectrum extends in some cases up 
to the infrared range (IChatv et al. Il2003t iKalemci et al. I 
[2001 iMighari et al. II2007D . At higher energies, from the 
optical to the soft X-ray frequencies, the emission is dom- 
inated by a thermal irradiated disk and the companion 
star, while at even higher energies, hard X-ray to soft 
7-rays, a power law spectral distribution is typically ob- 
served, with power law index 0.5 S, 5-, 1-1 and var iable 
exponential cutoff (jMcClintock &: Remillard I l2006l and 
references therein). 

Two classes of models are often used in interpreting 
the non-thermal part of the X-ray s pectru m : accretion- 
based models (e.g., iTitarchukI 1X9941: lEsin et I 
19971 iPoutanen I 119981: iCoppi I 119991: iZdziarskil 120001 : 



emission (see, e.g.. iMirabel fc Rodriguezl 119991 : [Fender I 
I2OO6I ). This idea had gained support by the correla- 
tion found bet ween the observed radio and X-ray emis- 
sionjgannikaineiL^Lay 119981; ICorbel et ani2000L 120031 
iGallo e it a.l. \\2(M\ . 

A model that explains the flat radio spectra observed 
by the VLBI in several compact o bjects as due to emis- 



sion from jets was suggested by Blandford & KoniglJ 
(1979). In this seminal work, |Blandford & KonigT] 
19791) showed that a flat radio spectrum can naturally 



.Esin et al. I 200l|) an d alternatively jet-bas ed models 
iMarkoff et al. " 



120011 [2003 20_05; Bosch-Ramon et aTl 

120061 iParedes et al. I [200I " iKaiser 1120061 )7" Indeed, as 
jets are widely believed to be the source of the observed 
non-thermal radio spectrum, it is possible that they also 
play a signiflcant role as a source of the hard X-ray 
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be obtained due to the change of the plasma conditions 
along the jet (in particular, the decay of the magnetic 
field and the decrease in the particle number density). 
This follows from the self absorption of synchrotron pho- 
tons, which produces a pronounced spectral break at 
vthick, the frequency below which the optical depth be- 
comes larger than unity. The flat radio spectrum follows 
the dependence of this frequency on the changing plasma 
conditions along the jet. 

In spite of their successes in reproducing the observed 
flat radio spectra, jet models are sti ll incomplete. In 
some of the models (e.g. , Blandfo rd fc Konigl I 1197^ ; 
IMarkoff et d^ |2001| . 120051 : IBosch-Ramon et al. I 12006 ) 
adiabatic, as well as radiative energy losses, are as- 
sumed to be fully replenished by an unknown, continu- 
ous re-acceleration process along the entire jet. While 
re-acceleration can result from in t ernal shocks within 
the jet (e.g iKaiser et al. I l2000t ISpada et al \ 120011 : 
iJamil et a/. 1120081 ). there is no a-priory reason to as- 
sume that it replenishes completely both the adiabatic 
and radiative energy losses, since the physical origin 
of these phenomena are d ifferent. A model developed 
by iHiellming fc Johnston I (jl98S) considered only adia- 
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batic, but not radiative energy losses of the electrons. 
En ergy losses of relativistic ele c trons was considered 
by iGeorganopoulos fc Marscher I (|1998l ). However, this 
model considered only the optically thin part of the spec- 
trum, and electrons propagation close to the line of sight. 
A comprehensive study of jet emission which in cludes en- 
ergy losses was carried out bv lRevnolds I ()1982D . although 
in this work the effect of self absorption on the electrons 

energy spectrum was not considered^ 

A different approach was used by Kaiser ' (2006) , who 
suggested a model in which only a single acceleration 
episode takes place, and radiation and a diabatic energy 
losses are not replenished. In this work, iKaiser I (|2006l ) 
pointed out that due to the decay of the magnetic field 
along the jet and because of self absorption effects, flat 
radio spectra can naturally be obtained. ,Kaiser , (,2006,) 
considered both conical jet in the ballistic case (i.e., ne- 
glecting adiabatic energy losses) , and a scenario in which 
adiabatic energy losses are considered. In both cases flat 
radio spectrum could be obtained, albeit with specific 
jet geometry (non-conical jet is required when adiabatic 
energy losses are considered). However, fitting the spec- 
trum of Cyg X-1 in this model could only be done with 
a very narrow jet, both in the ballistic and adiabatic sce- 
narios. 



In the model bv lKaiserl (j2006[ ) presented above (as well 
as some of the former models), a power law distribution 
of the accelerated electrons in the entire energy range 
was assumed. The break predicted in the spectrum is 
thus only due to the change from optically thick plasma 
at low energies, to optically thin above vthick- While the 
details of particles acceleration inside jets are uncertain, 
models of particles acceleration in shock waves predict 
that the power law distribution of electrons exist only at 
high energies, while a significant fraction, perhaps most 
of the electrons maintain lo w energy Maxwellian (ther- 
mal) distribution (Blandfo rd k Eichler I Il987t lAxford I 



SpitkovskvlHooil . Thus, 



119941: lAmato fc Blasi 200 
an inherent peak in the electrons energy distribution is 
expected, at typical Lorentz factor 7„iin- As a result, 
a peak in the synchrotron spectrum is expected at fre- 
quency Vpeak , due to cmission from electrons at the peak 
of the Maxwellian distribution. If Vpeak < vthick then this 
peak is obscured; however, if Vpeak > vthick, which, as we 
show below, is the case for plasma parameters that are 
plausible for jets in black hole binaries, then two breaks 
in the spectrum are expected, at vthick and at fpeafe- 

In addition to the characteristic low energy cutoff in 
the electrons distribution (attributed to electrons at the 
peak of the Maxwellian) any physical acceleration mech- 
anism accelerates electrons only up to a maximal Lorentz 
factor, denoted here by 7max- This introduces a maximal 
frequency t'max u p to which synchrotro n photons can be 
emitted (see, e.g.. iMarkoff et al. Il200l[ l. 

The uncertainty that exists in the nature of the ac- 
celeration process is reflected in an inherent uncertainty 
in the distribution of the energetic electrons. A power 
law energetic tail is often assumed; however, in addi- 
tion to the uncertainty in the fraction of particles that 
are accelerated to the energetic tail, there is an un- 
certainty in the spectral index of the power-law accel- 
erate d particles. Fitting d ata in supernovae remnants 
(e.g. iLazendic et al. I [20041) and gamma-ray bursts af- 
terglow (jGalama et al. 1119981 IWiiers fc Galama 1 119991 : 



iFreedman fc Waxman I l2001h reveal power law distri- 
bution with a nearly universal index, p — 2.2 ± 
0.2, which is consistent w ith theoretical expectations 
()Blandford fc Eichler |[T987[) . The X-ray spectra of many 
BHXRBs in the low/ha rd state shows photon inde x that 
average at < T >~ 1.7 (jMcChntock fc Remillard1l2006 ). 
If interpreted in the framework of the synchrotron emis- 
sion model, this implies electron spectral index p ~ 2.4, 
and in several cases even higher. The power law indices 
inferred from fitting synchrotron emission to the X-ray 
spectra of BHXRBs are thus consistently higher than 
those obtained in other objects. 

An additional inherent uncertainty in the study of 
emission from BHXRB's is related to the production and 
structure o f the jet. Th e jet m ay be conical, as assumed 
in Blandf ord fc KonigTI (|1979l ). or confined. Th e jet ge- 
omet ry can strongly affect the radio spectra (jKaiser I 
l2006f ). Another uncertainty exists in determining the 
microphysical processes in the jetted plasma, which gov- 
erns the evolution of the magnetic field and the proper- 
ties of particle acceleration to high energies, either via 
shock waves or other mechanisms (e.g., magnetic recon- 
nection). The energy in the magnetic field is sometimes 
assumed to be in equipartitiq n with the energy carried by 
the electrons (jLongair 1119941 ). However, several authors 
sugge sted the hypothesis of P o ynting-flux dominated jets 
(e.g.. lUstvugova^F "aTl 120001 : [Giannios fc Spruit I l2006t 



Komissarov et al. Indeed, in many fits (e.g., 

Markoff et al. 1 12005: Gallo et al. I l2007t iMighari et aU 



2007f ). the energy in the magnetic field is kept as a free 



parameter, and the fitted value is somewhat higher than 
equipartition value. As we show below, the value of the 
magnetic field affects the entire shape of the spectra in a 
non trivial way. 

In this paper, we combine the idea of a single acceler- 
ation episode first proposed by Kaiser (2006), with the 
recent results on the distribution of accelerated particles. 
We study here a model in which particles are acceler- 
ated only once at the base of the jet, and cool via syn- 
chrotron emission and possible adiabatic energy losses as 
they propagate along the jet. The magnetic field decays 
along the jet from its initial value at the jet base, Bq, 
which is taken here as a free parameter. We assume that 
the acceleration process produces a low energy cutoff in 
the energetic electrons energy distribution, in the form 
of a low energy Maxwellian, accompanied by a possible 
energetic tail. We show that the inclusion of this inher- 
ent cutoff leads to a variety of complex spectras, whic h 
are very different than that obtained bv iKaiserl ()2006f ). 

This paper is organized as follows. We first introduce 
the basic assumptions of our model in fj2l We solve the 
equations that describe the cooling rate of the energetic 
electrons in H2.3[ These equations hold the key to the rest 
of the analysis. We discuss in ^Slthe general properties of 
the resulting spectra, and show that qualitatively differ- 
ent results are obtained for wide and narrow jets (which 
are defined there) . We then study the resulting spectra in 
the different scenarios in details. In ij4]we study the basic 
model of synchrotron emission (neglecting adiabatic en- 
ergy losses) from a Maxwellian distribution of electrons. 
We extend our study in ^to include power law distribu- 
tion of energetic electrons. In |j6]and Sj7]we consider the 
effect of adiabatic energy losses on the obtained spectra 
for initial Maxwellian and power law distribution of the 
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energetic electrons, respectively. We then consider in ij8] 
the scenario of narrow jets. We summarize and conclude 
in iJHl emphasising the key properties of our model in 
view of the existing broad-band data sets. Details of the 
numerical model that is used in producing the spectra in 
Figg. [2]-[Il]are given in appendix lAl 

2. THE MODEL: BASIC PHYSICAL ASSUMPTIONS 

Our basic je t geometry is similar to the one treated by 
iKaiser I (|2006[ ). We consider a jet centered around the 
X axis, with rotational symmetry around that axis. De- 
noting by xq the position at the base of the jet, the jet 
radius with respect to the x-axis at x > is a func- 
tion of X only, r = r(x). In the following, we assume a 
parametric dependence r{x) = ro{x/xo)"'^'K Here, tq is 
the jet radius with respect to the x-axis at xq, and a jet 
defines the jet geometry: ajet = 1 corresponds to a con- 
ical jet. Since the exact shape of BHXRBs jets are not 
precisely known, we keep ajet as a free parameter. In the 
calculation below, we assume that the observer is located 
at an angle 6°^ ^ to the jet axis. This assumption al- 
lows us to simplify the radiative transfer calculations of 
synchrotron photons along the jet. 

We assume that the bulk motion Lorentz factor of the 
plasma inside the jet is constant, and equals 7j. This 
assumption can only be justifi ed for a ballistic jet ( a,et = 
1). However, it was shown bv lBlandford fc Rees I (|1974f ) 
that the Lorentz factor of the bulk motion of a flow in a 
jet is proportional to px where is the pressure of 
the external medium that confines the jet. Assuming a 
shallow pressure gradient, the Lorentz factor of the flow 
is not expected to vary significantly along the jet. A 
similar argument holds for magnetically confined jet. 

The plasma in the jet originates from the accreting 
disk, hence its density depends on the disk accretion 
rate , Mdisk- Following iFalcke fc Biermann I (|1995D (see 
also lMarkoff et al. 11200 ll ). we assume a fraction qj of the 
accreting matter to be injected into the jet. These as- 
sumptions imply that the (comoving) energy density at 
the base of the jet is given by uq — qj MdiskC^ / t^Tq^j f3jC, 
where /3j = (1 — 7^"^)"'^^^ is the bulk motion velocity, and 
ro ^ xo assumed. 

We consider here an acceleration process that occurs 
entirely at xq . This assumption is of course too simplis- 
tic. Nonetheless, we use it here in order to demonstrate 
the variety of spectra that can be obtained, and the spec- 
tral dependence on the uncertain value of the magnetic 
field, the jet geometry and the distribution of energetic 
particles. Given the existing uncertainty in the nature of 
the acceleration process, we use two models for the accel- 
erated electrons energy distribution, which we consider 
as the two extreme cases. 

The first is a (relativistic) Maxwellian distribution, 
nei{'y)d'y = AP'y'^ exp{—^ /9ei)d'y, where 7 is the Lorentz 
factor associated with the electrons random motion (not 
to be confused with the Lorentz factor 7^ of the bulk 
motion jet flow), /? = (1 — 7~^)^^^ is their random veloc- 
ity and A is a normalizaton constant whos exact value is 
determined below. The characteristic, normalized elec- 
trons temperature, O^i at the acceleration site x = Xq is 
determined as follows: we assume that the plasma in the 
jet contains both electrons and protons. The accelera- 
tion process (e.g., a shock wave) dissipates the kinetic 



energy of the flow, and redistributes the proton energy 
between the electrons and the protons. Hence, the elec- 
trons get some fraction Ce of the equipartition value of 
the total (electrons + protons) energy. Therefore, the 
average Lorentz factor associated with the electrons ran- 
dom motion is 7min = ee(TOp/™e)- Since in a Maxwellian 
energy distribution the average Lorentz factor is related 
to the temperature via 7min = 36'ez , we conclude that the 
characteristic electrons temperature at the base of the 
jet is 

/ 6 \ 7TI 

9ei {x^xo) = Oei.o = hr — = 612ee.o. (1) 

V 3 / TOe 

We adopt here and below the standard convention Q = 
lO^Qx in CGS units. 

Since a large number of pairs is not expected to be 
created, the total number density of electrons is similar 
to the protons number density in the plasma. At the 
base of the jet, it is equal to Uei, tot ~ uo/mpC^. The 
normalization of the accelerated electrons distribution is 
given hy A — ne;^ tot/26'g; q' where the factor 2 comes 
from integration over the electrons energy distribution: 
/r nei(n)d'y = A0ei,ok2iOeLo)^ where k2iO~^o) ~ 26*^; q is 
Bessel k- function of second order, and the last equality 
holds for 9eLo » 1. 

The second electrons energy distribution that is consid- 
ered here is a power law distribution above 7min, while 
maintaining a Maxwellian distribution of the electrons 
population at lower energies. This is motivated by the 
theoretical works discussed in fJT] In this scenario, the 
power law index of the accelerated electrons energy dis- 
tribution p = — (ilognei(7)/dlog(7) is taken as a free 
parameter. The electrons are accelerated up to maxi- 
mal Lorentz factor 7niax, whose value is determined by 
equating the acceleration time and the synchrotron loss 
time (see discussion in §2.1|) . An exponential cutoff in 
the electrons energy distribution above 7max is assumed. 

The two distributions considered here rely on strong 
theoretical basis and are consistent with modeling spec- 
tra from jets in gamrn a-ray bursts (GRBs; see, e.g., 
iPe'er fc Waxman II200I ). We note though that alterna- 
tive acceleration models ex ist. E.g., in modelin g emis - 
sion from jets in Blazars, ICelotti fc Ghisellini I (j2008l ) 
obtained good fits to the data with electron distribu- 
tion that is consistent with power law at low energies, 

7inin > 1- 

The magnetic field evolves along the jet. We consider 
here a scenario in which the magnetic field is dominated 
by the toroidal component, B = Bcj, oc r~^. Therefore, 
B = Boir/ro)"^ = Bo{x/xo)~'''^* ■ The value of the mag- 
netic field at the base of the jet, Bq is a free parameter. 
One way to quantify it is by assuming that the magnetic 
field carries some fraction es of the dissipated kinetic 
energy: ub — ^bUq, resulting in Bq — (STresUo)^^^- An 
equipartition value of the magnetic field is therefore ob- 
tained by setting = 1- 

As the electrons propagate along the jet at x > xq, 
they cool via synchrotron emission and possible adiabatic 
energy losses. In this work, we do not consider other 
physical processes (such as, e.g., Compton scattering, or 
production of pairs). This is due to two reasons: first, for 
plausible parameters describing jets in BHXRB s, other 
phen omena are much less significant (see, e.g. iKaiser I 
I2OOI . Second, the aim of this paper is to show that syn- 
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chrotron emission combined with adiabatic energy losses 
can lead by themselves to a very large variety of possi- 
ble spectra. When adiabatic energy losses are included, 
we use the precise formula, which is correct in the ultra- 
relativistic as well as the non-relativistic regimes. 



jet is 



91og(7/3) 
dt 



i aiog(AT/) 
'3 dt 



(2) 



which leads to 7/3 oc {AV)~^^^. Here, AV is the volume 
element occupied by the particles; since the particles are 
assumed to propagate at constant velocity in the x di- 
rection, AV oc r(x)^. We thus find that adiabatic en- 
ergy losses result in 7/? oc r^^/"^. For energetic electrons, 
7^1, this formula asymptotes to the more familiar 
form, 7 oc 7-^2/3^ 

In calculating the observed synchrotron flux, one needs 
to solve the radiative transfer equation along the line of 
sight. We use the standard assumption that the radiation 
field is isotropic in the comoving frame of the fiuid. In 
this work we limit ourselves to consider observer location 
at high angle to the jet propagation axis. This enables 
us to carry these calculations analytically. The results 
of the analytical approximations are then checked with 
the more precise numerical calculations, whose details 
are described in appendix i jAl 

We calculate the emission in the direction perpendicu- 
lar to the X axis by splitting the jet into small segments 
of length dx, and carrying the radiative transfer calcu- 
lations at each segment independently. The surface area 
of a segment is 2'Kr[x)dx. The calculation can thus be 
considered as a first order approximation, although the 
errors are not expected to be large as long as the angle 
between the jet axis and the observer, 9°^ ^ 0. In the 
calculation of the observed intensity below, we further 
omit the Doppler factors for the approaching and reced- 
ing jets, = [7i(l T Pj cos6'°'')]~^. For arbitrary angle 
to the line of sight 9°^, the Doppler factors arc within the 
range < 6^ < (27,)-! (5^^" - 7r/2] =77!), and 
are therefore not much different than unity for mildly 
relativistic jets, 7j 1- 

2.1. Characteristic frequencies of synchrotron emission 
from electrons at the base of the jet 

2.1.1. Maxwellian distribution of the accelerated electrons 

We first consider a relativistic Maxwellian energy dis- 
tribution of electrons at the base of the jet, nei{j)d'y = 
y4/37^ exp(— 7/0e/.o), with a normalized temperature 
given by Eq. [TJ In calculating the characteristic 
break frequencies of the observed emission in our model, 
we use the same par ameters values that were used by 
iMarkoff et al. I (|2001l ) in fitting the broad band spectrum 
of XTE J1118+480. This gives possible values of the free 
model parameters that can be used for illustrative pur- 
poses. 

We therefore assume a disk accretion rate Mdisk = 
3 X 10"* Mq yr"i around a Mbh = 6M0 black hole, an 
efficiency in matter injection into the jet qj = 10"^ and 
bulk motion Lorentz factor of the matter inside the jet 
7j — 2. We take the jet base to be at distance xq = 45rs 
from the black hole, and the jet radius at its base to be 
ro = lOr^, where r^ is the Schwarzschild radius. With 
these assumptions, the energy density at the base of the 



Uo 



qjMdiskC 10' 



,26 



—uq erg cm 



(3) 



where uq = gj,_2(Mdisk/lO'''^M0 yr~-^)(7j/2) is a dimen- 
sionless quantity, whose value depends on the fiow pa- 
rameters. For the nominal values taken here, the nor- 
malization of the electrons distribution, Aq = A{x — 
xo) = 5 X 10^ -So T-Q^J e^Q cm"'^, and the magnetic field 

at the base of the jet, i?o = 9 x 10* Uq^^ r^l e]/^3 G are 
readily determined. Here, ^0,1= ro/lOrj. 

The peak frequency of the spectrum is determined by 
synchrotron emission from electrons at the peak of the 
Maxwellian. At the base of the jet. 



^peak^O 



l^peakix ^ Xo) = 
1.4 X 10" U^' 



3^ei.o qBp 
47r mgC 
-1.2 1/2 TT 
0,1 ^e,0 ^B.-3 



(4) 



Synchrotron emissivity from a Max wellian distribu- 
tion of electrons was calculated by I Jones fc Hardee 1 
(fT979f >. For z — /2vp(,ak ^ 1, it was found that 
< 1) = (4/9)(^6'3,g3So/mec2)(z/2)i/3, while for 
z 3> 1 the fiux decays exponentially, iu{z 3> 1) = 
(7r/4)(Aeg3Soz/m,c2)exp[-(3/2)(2z)i/3]. 

For a propagation direction perpendicular to the x 
axis, the photons travel distance r{x) before escaping 
the jet. The optical depth to sync hrotron self absorption 
of th ese photons is given by (e.g.. iRvbicki fc: Lightmanl 
[T979h 

= (5) 
Zv^mf.c'^9ei 

The break frequency is defined as the frequency below 
which the optical depth is larger than unity, vtMck = 
v\r^=i- Taking vthick < Vpeak (i.e., z < 1) at the jet 
base, using Eqs. [4]and[5l one finds that at the base of 
the jet 



Vthick = Vthick{x — Xq) 



9 m^c^ 
9^ m^m^c^ 



3/5 



1/5 



1/5 



1 

ci.O 



(6) 



2.1x1013 ^^^/^o:leeJ4(-3Hz 



The condition Vpeakfi > vthick,o can thus be written as a 
requirement on the energy carried by the electrons, 



^ c o w in-2 -1/10 -1/10 



(7) 



The weak dependence on the parameters of the flow and 
on the magnetic field, implies that for a wide range of 
parameters, the right hand side of Eq. [7] is smaller than 
unity. Thus, for near equipartition value of the energy 
carried by the electrons tefi ^ 1, the peak emission fre- 
quency at the jet base is above the break frequency. 

2.1.2. Power law distribution of the accelerated electrons 

Our second acceleration scenario considers power law 
distribution of the energetic electrons above the peak of 
the MaxweUian: 726/(7)^7 oc 7"^ for 7inin < 7 < 7max- 
The power law index p of the electrons energy distri- 
bution is taken as a free parameter. We estimate the 
maximum Lorentz factor of the accelerated electrons by 



Jets in XRBs 



5 



equating the acceleration time, tacc — E/cqBo^ with 
the synchrotron cooling time, Qmlc' / Aq'^ Bq'^ , to ob- 
tain 7,„ax = (3/2)r7ieC^/('Z'^^o)^^^- Synchrotron emission 
from these electrons is expected at 



^max,0 



3 iL^qBo 27 mec3 



Att rrieC 



IGtt 



5.7 X 10- Hz. (8) 



Note that this value is independent on the strength of 
the magnetic field at the acceleration site. 

The minimum Lorentz factor of the power law distribu- 
tion, 7niin is determined self consistently once the power 
law index p, the maximum Lorentz factor 7max and the 
total number and energy densities of the accelerated elec- 
trons are known (see further explanation on the numeri- 
cal model in fJX|. Below 7min, a Maxwellian distribution 
is assumed, with temperature 0ei,o — 7min/2, chosen such 
that the low energy Maxwellian smoothly connects to the 
power law distribution at higher energies. The inclusion 
of a low energy Maxwellian distribution implies that the 
results derived above in Eqs. S]- [7] hold in this scenario 
as well. 

2.2. Observed flux from a jet segment 

For an observer located at high angle to the a;-axis, a jet 
segment has a surface area 2TTr{x)dx. The flux observed 
from a jet segment of length dx is thus 



dF„ = 



■ix)^j^ 

2d2 



(1 



')dx, 



(9) 



where d is the distance, and we omitted the dependence 
on the Doppler factors. 

For a Maxwellian distribution of electrons, the ob- 
served spectrum from a jet segment can therefore be de- 
scribed by a broken power law with two characteristic 
frequencies: for v <^ Vthicki ''"i^ ^ Ij and 



2d2 n 



dx — 



d2 



nieOeii'^dx, (10) 



is proportional to i/^. For vthick ^ ^ ^ i^peak, the optical 
depth T,y < 1 and 



jijdx 

It 



1/3 



AFu{vthick I^V < I^peak) 

. . . ^^^^ 

is proportional to ly^^'^. For v > i^peak, the emissivity ji, 
decays exponentially, and therefore an exponential cutoff 
in the observed flux is expected. 

In this work we focus on a scenario in which Cg is close 
to equipartition. However, we note that if the electrons 
carry a much smaller fraction of the energy such that the 
condition in Eq. [7] is not fulfilled, then vthickfl > fpeafc,o- 
Moreover, even if eg — 1, as we will show below both 
i^peak and Vthick Can vary along the jet in different ways. 
As a result, there exists a transition radius Xtrans above 

which Vth%ck{x > Xtrans) > l^peak{x > Xtrans)- OuCC 

this occurs, the emission from a jet segment is composed 
of a thick part below vthick where dF^ oc v'^, and an 
exponential cutoff at higher frequencies. 

When power law distribution of electrons is considered, 
in order to calculate the emissivity and flux one needs to 
specify the fraction of electrons that are accelerated to 



the high energy tail above the Maxwellian. For demon- 
stration purposes, we consider here complete acceleration 
to the energetic tail, while maintaining the low energy 
Maxwellian and high energy cutoff in the distribution. 
Thus, the power law distribution smoothly connects to 
the peak of the Maxwellian. In this case, calculation 
of the frequency dependence of the flux in Eqs. [H - [TT] 
hold (A slight modification though exist in the overall 
normalization of the flux). 

We can thus conclude that there are two main differ- 
ences in the emissivity from a jet segment in a power law 
scenario and the pure Maxwellian distribution considered 
abov e. The first difference is tha t standard calculations 
(e.g. jRvbicki &: Lightman]|1979D show that above Vpeak 
(and below t'max)7 there is a power law decrease of the 
flux, with power law index {p—l)/2. An exponential cut- 
off exists only above fmax, while in a pure Maxwellian 
distribution an exponential cutoff exists already above 

Vpeak • 

The second difference is relevant if vthick{x) > 
Vpeak{.x). For powcr law distribution of electrons, the 
thick part of the spectrum (y < vtMck) is characterized 
by F^ (X fe.g-. lRvbicki fc LightmMnil979D . while a 
power law F, oc i/(P-i)/2 is kept above vthick- Here, how- 
ever, we expect a somewhat more complex spectral shape 
at low frequencies i/ < lytMck- for v < Uthick, F^ cx v^/"^, 
while at much lower frequencies, the flux is Fjj cx due 
to the inclusion of a low energy Maxwellian. 

2.3. Electrons energy loss along the jet 

As the electrons propagate along the jet, they cool by 
synchrotron emission, as well as possible adiabatic energy 
losses. As a result, the characteristic frequencies derived 
above for x = xq vary along the jet. We calculate in this 
section the change in the electrons energy distribution 
due to synchrotron and adiabatic cooling. 

As they propagate along the jet, the electrons cool 
via synchrotron emission, and their momentum decreases 
at a rate d(7/3)/c?i|sy„c = ~{A/i)caT I {■meC^)uB{ll3)'^ ■ 
Here, ax is Thomson's cross section and ub — i?^/87r 
is the energy density in the magnetic field. Along the 
jet, the magnetic field decays as B{x) — Bo{x/xo)~°-^'* . 
Since the bulk velocity of the electrons propagation is 
taken to be constant, the electrons position along the jet 
(x) axis is related to the (comoving) time they spend in 
the jet by a; = jjPjCt; thus, the electrons reach the accel- 
eration site (at the base of the jet) at time to — xo/^jPjC. 

Adiabatic energy losses (see Eq. [2]) are quantified by 
d-i"^ (3) / dt\^^ = —{2/3){ajet/t)jP. Considering both syn- 
chrotron and adiabatic energy losses, the electrons mo- 
mentum dec ay (in the bulk motion comoving frame) is 
governed by (jLongair 111994; Kaiser 2006) 

di^^_AaT_Bl/iy^''^ 
3meC 87r \to/ 

(12) 

The solution to Eq. [T^] was derived by iKaiser I ()200i 



dt 



7/3(i) 



7/3(to) X (i/to)-^'^-'/^ 



1 



+ ,.,nX.::/3-l, lf^it0) [1 - (tAo)l-«'^-'/^] 

(13) 

If adiabatic energy losses are neglected, the second term 
on the right hand side of Eq. [T2l drops, and the solution 
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-^l(3{to) [1 - (tAo)i 



-2a,- 



(14) 

Equation [14] reveals an important result, which was 
first derived by Kaiser (2006): if adiabatic energy losses 
are not included, for a^et > 1/2, the Lorentz factor of 
the electrons asymptotes at t ^ oo to a value larger 
than unity. This implies that some part of the electrons 
initial energy is maintained and not radiated. The ori- 
gin of this counter-intuitive result is the decay of the 
magnetic field along the jet. It provides a natural, pos- 
sible mechanism for maintaining the electrons (asymp- 
tote) temperature, a s is req uired by, e.g., the model of 
iBlandford fc Konigl I ()1979D . without the need for par- 
ticle re-acceleration along the jet. This fact also al- 
lows us to define wide jets in this case as jets for which 
Ojet > 1/2. 

We present in Fig. [1] several examples of the electrons 
momentum temporal decay due to synchrotron and adia- 
batic energy losses (Eq. I13p and pure synchrotron cooling 
(Eq. [T4|) . In producing the results, we consider a conical 



1), and four different values of the magnetic 



jet (fljet 

field. Bo = 10*, 104■^ 10^ 10^-^ G. An additional plot il- 
lustrates the scenario of a narrow jet (ojej = 0.1) and 
Bq = 10^-^ G. Asymptotic behavior of the electrons mo- 
mentum at late times, t ^ to, is clearly seen (and easily 
derived from Eqs. [T3 lfT4|) . When adiabatic energy losses 
are considered, for wide jets, aj^t > 3/8 the electrons 
momentum decays asymptotically as 7/3 oc f-'^o-j^t/^ ^ Pqj. 
pure synchrotron cooling and aj^t > 1/2, the electrons 
momentum reaches a constant value. Interestingly, we 
find from Eqs. [T51 and [HI that in narrow jets, ajet < 3/8, 
the electrons asymptotic decay law is similar in both sce- 
narios, i.e., independent on the inclusion of adiabatic 
energy losses. In this case, for t to, both scenar- 
ios result in electrons momentum asymptotic decay law 

For strong value of the magnetic field, the asymptotic 
behavior of the electrons momentum is reached after a 
very rapid decay of the electrons initial energy, or mo- 
mentum. This rapid decay results from an extensive syn- 
chrotron emission in the strong magnetic field close to 
the jet base. The results presented in Fig. [T] illustrate 
an important point that will be extensively used in the 
analysis in the following sections: in a strong magnetic 
field, the initial rapid decay in the electrons momentum 
takes place on a very short time scale, t/to < few, which 
is translated to very short spatial scale - the rapid cool- 
ing (and most of the radiation) occurs very close to the 
jet base. 

This enables us to obtain an analytical approximation 
to the observed flux by splitting the jet into two sepa- 
rated regimes. We first calculate the emissivity during 
the rapid cooling phase, which takes place close to the jet 
base, hence the magnetic field can be approximated as 
constant during this phase. At a second step, we calcu- 
late the emissivity during the rest of the electrons prop- 
agation along the jet, assuming that the electrons tem- 
perature evolution follows its asymptotic behavior. The 
approximate analytical results are compared with the ex- 
act results obtained numerically. 
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Fig. 1. — Temporal decay of the electrons momentum, 7/3(t) as 
a f unct ion of time. Shown are adiabatic -|- sync hro tron (solid lines, 
eq. I13I I and pure synchrotron (dashed lines, eq. I14I I scenarios. Four 
values of magnetic field are considered (from top to button) , Bq = 
10*, 10'*-5,105, lO^-^ G. The dotted lines (light blue) show the decay 
of the electrons energy for narrow jet, ajet = 0.1 and Bq = 10*'^ G. 
In all the other cases we consider conical jet (ojet = 1)- The values 
of the flo w para meters used in determining to are the same as those 
used in i|2.1.1l resulting in io = x^j'^jfijC = 1.5 X 10~^ s, and 
initial electrons momentum -fl3^i{to) = mp/nie = 1836. For pure 
synchrotron and ajgt > 1/2, the momentum asymptotes at t 2> to 
to a finite value, while for synchrotron and adiabatic losses, at 
t ^ to, 7/9(t) c< t~^°-3i!t/^. An important result is that the initial 
rapid decay occurs on time scale t/to ^ f^w. See text for details. 

3. GENERAL PROPERTIES OF THE SPECTRA 

As we will show in iJH - ^[51 a wealth of spectra can be 
obtained, whos details depend on the various model as- 
sumptions. In particular, we will show that the observed 
spectra is very sensitive to the strength of the magnetic 
field. This dependence is found to be highly non-trivial. 
In addition, both the jet geometry and the initial par- 
ticle distribution significantly affect the obtained spec- 
trum. Therefore, before discussing the various possibili- 
ties in details, we first discuss in this section some general 
properties of the emission along the jet. Here we focus 
on global properties of the observed spectra, which are 
direct consequences of our models assumptions as stated 
in 321 The detailed spectra expected in the various sce- 
narios will be discussed in the following sections. 

3.1. Critical values of the magnetic field at the base of 

the jet 

We showed in H2.'3\ that for strong magnetic field at 
the jet base, the electrons temporal behavior can be sep- 
arated roughly into two distinctive regimes (see figure [T]). 
At the first stage, the electrons rapidly cool due to syn- 
chrotron emission. This stage lasts a very short duration, 
of the order t/to ^5 f^w. Following the rapid energy loss, 
at a second stage the electrons cooling rate asymptotes 
to its terminal value. If only synchrotron emission is 
considered and a jet > 1/2, then the electrons maintain 
their energy, and their Lorentz factor, or momentum, 
becomes time independent, 7/3(t) oc t*^ at t ^ to. If adi- 
abatic energy losses are considered, then for a jet > 3/8 
the electrons momentum decays at late times t ^ to in 
accordance to 7/3(t) oc ^-^ajet/s^ -pov lower values of ajet, 
both scenarios yield a similar result, 7/3(t) oc p'^j't-^. 

The fact that the rapid decay phase exists only very 
close to the jet base, in regions where the magnetic field is 
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approximately constant (and equals Bq), combined with 
the fact that the accelerated electrons distribution has 
characteristic energies, enables one to define critical val- 
ues for the strength of the magnetic field at the jet base. 
As we will show below, these critical values represent 
qualitative transitions in the resulting spectra. The crit- 
ical values are define as follows. 

If the acceleration process produces a Maxwellian en- 
ergy distribution, then the energetic electrons have a 
characteristic Lorentz factor, 7min(^o) = 30ei.o- Due to 
the exponential decay nature of this distribution at high 
energies, there are effectively no electrons with Lorentz 
factor 7 ;» 7min- The cooling rate depends both on the 
electrons energy, or momentum and on the strength of 
the magnetic field at the jet base, Bq (see Eqs. fT2l~ [T4)) . 
One can therefore conclude that for very weak magnetic 
field, the rapid decay phase is insignificant and the entire 
electrons distribution is effectively unchanged during this 
phase. This case is illustrated in Fig. [TJ for Bo = 10^ G. 
If, on the other hand, the magnetic field is stronger than 
a certain value, the initial electrons cooling is significant. 

This fact allows us to define a critical value of the mag- 
netic field at the jet base, as a value above which elec- 
trons at the peak of the Maxwellian distribution signifi- 
cantly lose their energy by synchrotron emission during 
the initial rapid cooling phase. For the pure synchrotron 
case (neglecting adiabatic energy losses) and wide jets 
{a jet > 1/2), this critical value is formally obtained by 
taking the limit t ^ to in Eq. [TH and equating the 
second term in the denominator with unity. 



Bcr.l — 



1.6 X I0'^{2ajet - 1)^/^ 
x.-y^(7,/2)V2(:,o/45r,)-V2 q, 

3 X 10-^ {2aj 



-I- jet 



1 



(15) 



While the results in Eq. [15] are formally derived for the 
pure synchrotron scenario in wide jets, one finds that the 
time scale of the initial, rapid cooling is similar when adi- 
abatic energy losses are included (see Eqs. [131 [HI and 
Fig. [1]) and in the case of narrow jets. Therefore, the 
results derived in Eq. [TSlhold also in theses cases (math- 
ematically, the only difference is to replace the factor 
2a jet — 1 with a similar factor 8a jet /3 — 1 when adiabatic 
energy losses are added). We can therefore summarize 
that for Bo < Bcr,i, the initial rapid cooling is insignifi- 
cant. 

Our second acceleration model considers power law en- 
ergy distribution of the accelerated electrons above the 
peak of the Maxwellian and below a maximum Lorentz 
factor 7max- Equations [T3l and [Ml that describe the de- 
cay of the electrons Lorentz factor due to synchrotron 
emission and adiabatic cooling are independent on the 
electrons distribution, and hence are valid in this sce- 
nario as well. 

The inclusion of a maximum Lorentz factor of the elec- 
trons energy distribution 7,nax implies that one can define 
an additional critical value of the magnetic field, Bcr,o in 
a similar way to the definition of Bcr.i- It is defined as 
the value of the magnetic field below which electrons at 
7inax do not cool significantly by synchrotron emission. 



(16) 



Repeating the analysis carried above shows that 

BerO = 270{2ajet - 1)^/^ 

x(7,/2)2/3(xo/45r,)-2/3G, 

es.cr.o = 10-« (2a,et - 1)^/^ 

x(7,/2)4/3(xo/45r,)-4/3^,-V2^^. 

We thus conclude that for very low values of the mag- 
netic field, Bq < Bcr,o, the distribution of the power law 
accelerated electrons in the entire energy range is unaf- 
fected by the initial rapid cooling phase. If the magnetic 
field is stronger, Bcrfi < Bq < Bcr,i, then the highest 
end of the electrons distribution is cooled rapidly close 
to the jet base, while for electrons at the low end of the 
distribution, with 7 > 7inin, the initial cooling is inef- 
ficient. As a result, the main effect of the initial rapid 
cooling is a shift of the highest energy tail in the elec- 
trons energy distribution to lower energies. This result 
affects the observed spectrum at high energies, and will 
be discussed in [J5l 

For higher values of the magnetic field, the initial rapid 
cooling of the energetic electrons is significant. In the 
pure synchrotron case and a jet > 1/2, following the ini- 
tial rapid decay, the electrons Lorentz factor asymptotes 
to a constant value. Since in this case at t to the sec- 
ond term in the denominators of both Eqs. ll3l and [T4l are 
larger than unity, one finds a simple analytical expression 
for the energetic electrons asymptotic momentum. 



(7/3)/ = -/Pit » to) 
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(17) 



The electrons Lorentz factor at the end of the rapid cool- 
ing phase is therefore independent on its initial Lorentz 
factor, and is inversely proportional to Bq (see Fig. [ij. 
Similar to the discussion that followed Eq. [TSl we point 
that while Eq. [T71 was derived for the pure synchrotron 
case, a similar equation holds when adiabatic losses are 
included. 

If the magnetic field is in the range Bcr.o < Bq < Ber,i 
then the rapid cooling affects only electrons at the high 
end of the power law energy distribution, while if the 
initial distribution is a Maxwellian, the entire electrons 
distribution is unchanged. For stronger value of the mag- 
netic field. Bo > Bcr,i, electrons at and above the peak of 
the Maxwellian rapidly cool. Since the electrons Lorentz 
factor at the end of the rapid cooling 7/ is independent 
on its initial Lorentz factor, then in both the Maxwellian 
and power law acceleration scenarios we expect a quasi- 
Maxwellian distribution to be formed at the end of the 
initial rapid cooling phase (see the numerical results in 
the figures below) . This quasi-Maxwellian distribution is 
characterized by temperature 9eij — 7//3. 

During the initial rapid decay phase, the magnetic field 
is nearly constant. However, the characteristic frequen- 
cies Vpeak and Vthick depend on the electrons tempera- 
ture, and therefore vary as the electrons cool (see Eqs. 
HllHl). While the peak frequency fpeak decreases when the 
electrons lose their energy, the break frequency vthick in- 
creases. Therefore, if the electrons cool to below a critical 
temperature, 

' ^ \ 1/3 / 2 \ 1/5 

4 \ ' / TT^uogro \ on -i/io -1/10 

-B.-3 ' il^i 



\9BompC 



= 30 u, 



then at the end of the rapid cooling phase, vtMckj = 
vthick{.x = xo;Oei = Oeij) > i^peakj ■ The Condition 
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that the electrons cool to temperature 9eij < Ocr can 
be phrased as a condition on the strength of the mag- 
netic field at the base of the jet. Using Eq. [T71 one finds 
that this happens for Bq > Bcr,2, where 

Bcr,2 = 8.0 X 10"^ (2ajet - 1)^/^ 

X (7,/2)5/9 (a;o/45r,)-5/9 u-'^' rl{^ G, 

eB,cr,2 = 8.0 X 10-4 (2ajet-l)'°/' 

X (7,/2)io/9 (xo/45r,)- W9 ^-"/s ,20/9^ 

' (19) 

We thus conclude that if the magnetic field is interme- 
diate, Bcr,i < Bo < Bcr,2, then at the end of the initial 
cooling phase, still the peak of the emission is in the opti- 
cally thin region, i^peakj > ^thickj- For higher values of 
the magnetic field, at the end of the rapid cooling phase 
the peak of the emission is obscured, vthickj > i^peakj- 

3.2. Asymptotic variations of the characteristic 
frequencies along the jet and the transition 
frequency 

The change in the electrons energy distribution and the 
decay of the magnetic field along the jet imply that the 
characteristic emission frequencies, vtMck and Vpeak vary 
along the jet. Eq nations [T51 and [HI give the decay laws of 
the electrons energy along the jet, and thus enable to cal- 
culate the variation of these frequencies. Although these 
equations show complex dependence of the electrons en- 
ergy in time, the discussions in iJ2.3l and iJ3.1l show that 
following an initial rapid decay phase, the electrons decay 
laws at times t ^to asymptote to simple functions, that 
can be described analytically. Therefore, simple analyti- 
cal description of the characteristic frequencies temporal 
dependence is obtained in this limit. 

If adiabatic energy losses are ignored, then for ajet > 
1/2 the electrons temperature asymptotes to a constant, 
finite value. Thus, we consider the asymptotic temporal 
dependence dei oc oc r°. In this case, the electrons 
number and energy densities drop as n^i^ totir) oc r^^, 
u{r) oc resulting in a decay of the proportionality 
constant A{r) cx r~'^. Since the magnetic field also de- 
cays as B{r) oc r~^, one finds from Eqs. [4]and[6]that the 
characteristic frequencies decay in a similar way, 

Vpeak oc VtMck « OC x^"^'"* . (20) 

When adiabatic energy losses are considered, for ajet > 
3/8, at t > to the electrons temperature decays as 
9ei{t) oc i-2ajet/3 ^ ^-2/3 rpj-^g elcctrons number den- 
sity drops as tot('') r~'^, and therefore the electrons 
energy density decays as u{r) oc r^^^^. As a result, the 
characteristic frequencies decay according to 



'peak 
Vthick 



BOl 



„-7/3 



oc X 

oc (^352^3)1/5^^-^1 ^ ^-11/15 ^ ^-2a,et/3^ 



— 7ajat /3 



. (21) 

The different a;-dependent of the two frequencies imply 
that when adiabatic energy losses are considered, even if 
at the base of the jet Vpeak.o > vthick,o, then far enough 
along the jet, at x > Xtrans, vthick{x) > Vpeak{x). In 
order to estimate the transition distance x trans j we dis- 
criminate between two cases: 

(I) Bo < BcrX- In weak magnetic field, the electrons 
do not undergo the initial rapid cooling phase, and one 



finds that Xtrans/Xo = {Vpeakfi/Vthickflf^^"'''", 



Xtr 



Xo 



(250) 



l/Oje 



-3/(16ajet) 15/(8ajet) 3/(16ajet) 
) ^e,0 ^B,-3 



(22) 
At 

X trans ^ Vthick(.Xtrans^ — Vpeak{xtrans\ and both are 

equal to 



Vtr 



,(So < BcrA) = 3.7x10 



11 ,-.15/16 -1 --L^/o H7 
"0 '0,1 ^e.O ^B,-3 

(23) 

Thus, a change in the spectrum at vtrans is expected. 

(II) Bcr,i < Bo'- In this case the electrons rapidly cool 
before their temperature asymptotes. Thus, at the end 
of the rapid cooling, the ratio between the peak and 
break frequencies is Vpeakj /vthick.j = Vpeak,o/vthick.,o x 
{Oeij/Oei,oY- Following the initial rapid cooling the 
asymptotic behaviour of the two frequencies follows a 
similar decay law to the one considered above in Eq. dTJ 
and therefore one finds that the transition from optically 
thick to optically thin emission occurs at frequency 

Vtrans{Bo > B,r,l) = 1-1 X lO^^ (Sfljet/S - l^-l^/S 



-19/8 1/16 
'b,-3 



-53/16 -23/4 



1) 

-19/8 



>^Uo''"ra;i' ' (7j/2) 
x(:Eo/45r,)i9/«e^^/_^3^Hz. 

(24) 

By definition, if Bo > Bcr,2 then the transition to the 
optically thick emission occurs already at the base of the 
jet. 

For narrow jets, a jet < 1/2, the electrons tempera- 



ture decays as Oei oc x 



and their energy density 



drops as u oc x ^. As a result, the peak frequency de- 
cays as Vpeak{x) OC x^°'^^*~'^ , whilc the break frequency as 
vthick oc a;(-^"j=*+2)/5^ ^^(j ^Yig^^ ^j^g transition fre- 
quency depends on the jet geometry in a complex way: 



Vtransiajet < 1/2) « 10^^ X (7 X 10^' 



-1/2-3^2/10 -1 

^"o ''0,1 

^,2+3A2 ,1/2+3^2/10 „ 



(25) 



where Ai(ajet) = (15ajet - 10)/(12 - 24ajet), and 
^2{0'jet) = (12 — 2Aajet){iO'jet — 2)/5. The transition 
frequency varies as a power in Ai{ajf,t)- For jet geom- 



etry a jet ^ 1/3, one finds that vtr 



10^3 Hz, how- 



ever Vtrans drops sharply for higher values of a jet'- for 



ajet — 0.40, Vtrans 
Vtrans ~ 10^ Hz. 



10^" Hz, while for a 



jet 



0.45, 



3.3. Observed break frequencies from synchrotron 
emission along the jet 

The analysis carried above implies that the observed 
flux from the entire jet (as opposed to a single jet seg- 
ment) is characterized by several break frequencies. We 
present here a general discussion on the nature of these 
break frequencies. This is done in order to establish a 
general basis for the detailed discussion on the various 
models (Maxwellian vs. power law, pure synchrotron vs. 
adiabatic, etc.) discussed in the following sections. 

We are able to identify five observed break frequencies. 
However, some of these frequencies are pronounced only 
under certain conditions. These frequencies are: 

(I) Vpeak,o (Eq- HI). This break frequency is inherent 
to any model that contains a low energy Maxwellian dis- 
tribution, and is therefore expected in all the scenarios 
considered. 
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(II) fmax.o (Eq. E]) is expected in models in which a 
power law distribution of the accelerated electrons ex- 
ist, since any acceleration model necessarily have a high 
energy cutoff. 

(III) When the magnetic field is higher than a mini- 
mum value, the electrons at the high end of the distri- 
bution rapidly cool. As a result, the cooled electrons 
at the base of the jet (with Lorentz factor 7/ given by 
Eq. I17p emit synchrotron radiation at characteristic fre- 
quency Vf = {3/4'!r){qBo/'mec){jl3)j: . In determining 
the observed frequency, one needs to discriminate be- 
tween the different cases. For a Maxwellian distribution 
and Bcr,i < Bq < i?cr,2, Vpeakj > Vthickj ■, and therefore 
emission at Vpeakj is not obscured. Thus, in this case 

y.e-J'Ll{^,/2f{x,/Abrs)-^\iz. 

(26) 

Note that while Eq. [221 was derived for a Maxwellian 
distribution, it holds for power law distribution as well 
as long as i?cr o < Bq < B^r i , when Q^i f is replaced by 
(7/3)/- 

For stronger magnetic field, Bq > Bcr,2, the electrons 
rapidly cool to temperature below the critical tempera- 
ture 9cr (see Eq. [T8| . Once this happens, the peak of 
the synchrotron emission is in the optically thick regime, 
and is therefore obscured: Vpeak.j < ^tuickj- The ob- 
served break in the flux therefore occurs at the transition 
frequency from the optically thin to the optically thick 
emission. By definition, this transition frequency occurs 
as the electrons cool to 6cr- Emission from electrons at 
this temperature peaks at frequency 



Vfast 

{Bo>B 

C 



39: 



l,cr- qBp 



4.0x1014 ul"\^\e%''^,Yi^. 



(27) 

(IV) When adiabatic energy losses are included, or for 
narrow jets a^et < 1/2, the analysis in i i3.2l shows that 
for Bq < Bcr,2 there exists a transition frequency vtrans, 
whos value was calculated in Eqs. [23] - [25l For high 
value of the magnetic field, Bq > Bcr,2, the emission 
peak is in the optically thick part of the spectrum already 
during the initial rapid cooling. As a result, in this case 
^trans — ^fast {Bq > Bcr,2), as is defined in Eq. [27l 

(V) The last transition frequency, which we denote 
here as viow is relevant in a scenario of power law en- 
ergy injection, when adiabatic energy losses (or narrow 
jets) are considered. Cooling of the electrons at the 
highest end of the distribution implies that the peak 
frequency of synchrotron emission from these electrons, 
Vmaxix) decays in a similar way to the decay of Vpeak- 
for wide jets with adiabatic energy losses, asymptotically 
t'max(a;) oc a;-7aj=t/3^ This decay is faster than the decay 
of Vthickj and therefore at a given distance along the jet 
xiow, these two frequencies become similar, and equal to 
viow Since 7max > 7min, uecessarily xiow > Xtrans, and 
hence viow < vtrans- We further note that if Bq > Bcr,i, 
then at the end of the rapid cooling, all the electrons 
with initial Lorentz factor above 7min cool to the same 
Lorentz factor, 7^. As a result, in this case viow — i^trans- 

We therefore concentrate in the regime Bq < Bcr,i- 
Equation [21] gives the decay law of vthick under the 



assumption of Maxwellian energy distribution, and is 
therefore valid as long as vtMck < i^peak, or equivalently 
X < Xtrans- At larger distances x > x trans, the decay 
law of vthick is modified and depends on the power law 
index p of the accelerated electrons. 

In calculating the decay law of vthir k in th is case, we 
largely follow the treatment bv iKaiserl (|2006f) . The elec- 
trons number density between 7min and 7max can be writ- 
ten as nei{'y)d'j = k{x)^~^d'y. Conservation of particles 
number along the jet can be written as k{x)"f~Pd'y = 

{AVo / AV)k{xo)^^^ d^o, where AV is the volume occu- 
pied by the electrons and the subscript '0' refer to the 
values at X = xq- Using Eq. [51 for adiabatic energy 
losses along the jet and using AV oc r^, one finds that 
k{x) = fc(xo)(r/ro)-(2p+4)/3. 

In this case of power law distribution, the emissiv- 
ity is > oc a nd the optical depth 

is oc ryfcB(P+2)/2iy-(p+4)/2 fe.g. JRvbicki fc LightmaiTI 
Il979| ). Using i^thick = v\r^=i, one finds 



Vthick OC 



2/(p+4) 



(28) 



where Ai(p) 



7p)/3{p + 4:). We emphasis that this 



decay law of I'thick is relevant only for x > Xtrans, while 



at smaller radii, the analysis carried in i^B. 21 holds. 

Repeating a similar analysis to the one carried in ^3.2\ 
one finds that for Bq < Bcrfli 



2.6 X 10^ 



-(3l2/16)(5-ii) -I5A2/W 



^0 



'04 



xe 



(A2/8)(15Ai-19) (A2/16)(2Ai+l) 



2.6 X 10^ ul^ 



.-2 ,-10 ,-0.4 
0,1 ^e,0 ^5,-3 



Hz, (29) 



Hz, 



where A2{p) = (7/3)/[(7/3) + Ai{p)], and the second 
line gives the approximate dependence for p in the range 
2.0 <p < 2.5. 
Similarly, for Bcr.o < Bq < Scr.i, 



10' 



r'r-f' (7,72)-^ (xo/45r,) 



2 ^1.5 



Hz. 

(30) 

Calculation of viow for narrow jets, a jet < 1/2 is 
straightforward, yet cumbersome. For 2.0 < p < 2.5 
and a jet ^ 0.1, we find i/iow ~ 10^^ Hz, while for wider 
jets, lyiow rapidly decays: i^iow ~ 10^° Hz for ajet = 0.2, 
and i/iow ~ 10 Hz for a jet = 1/3. As we will show in |j8l 
for narrow jets it is impossible to obtain a flat spectrum 
at the radio band, as suggested by observations. There- 
fore, we will give only brief explanation on the possible 
spectrum obtained in this scenario. 

We further discuss the role of these five transition fre- 
quencies on the observed spectrum in the following sec- 
tions, when we present a detailed analysis of the various 
spectras that can be obtained. The transition frequencies 
are marked in Figg. [2| - [9| when discussing the various 
possibilities. 

4. EMISSION FROM INITIAL MAXWELLIAN 
DISTRIBUTION OF ELECTRONS WITHOUT 
ADIABATIC COOLING 

In this section we present the results of a model that 
includes only synchrotron cooling of the electrons, i.e., 
we neglect adiabatic energy losses. While adiabatic en- 
ergy losses occur as the gas expands, both the formation 
and confinement mechanisms of jets are still not fully 
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understood. It could thus be that the expansion of par- 
ticles inside the jet is not purely adiabatic, and that the 
electrons are heated as they propagate along the jet (e.g., 
by lateral shock waves). Neglection of the adiabatic en- 
ergy losses can thus be considered as an extreme, yet still 
physically plausible scenario. 

We further concentrate in this section in a scenario in 
which the acceleration process results in a Maxwellian 
distribution of the energetic electrons. This assumption 
has both theoretical motivation, and is used for illus- 
trative purpose. From a theoretical perspective, as dis- 
cussed in fJT] the mechanism that accelerates particles to 
high energies (e.g., via shock waves) is not understood 
from first principles. While it is likely that some fraction 
of the electrons are accelerated to high energy power law 
tail above the peak of the Maxwellian, the fraction of 
electrons in the energetic tail is uncertain. Thus, we con- 
sider this scenario as being a physically plausible scenario 
of electrons acceleration process. In addition, studying 
this scenario is helpful in understanding the spectral de- 
pendence on the uncertain value of the magnetic field 
and possible jet geometries. We expand our model in ij5] 
to the cases of power law distribution of the energetic 
electrons and to consider adiabatic energy losses. 

Here and in sections ij5] - |J7] we focus on scenarios of 
wide jets. Narrow jet scenario is discussed in fJ5] 

4.1. Weak magnetic field: production of flat radio 
spectra in conical jets 

For weak magnetic field at the jet base, Bq < Bcr,i, 
and for wide jets, Ojet > 1/2, electrons cooling is in- 
significant in the entire jet (this scenario is illustrated 
in figure ^ Bq — 10* G). The insignificance of elec- 
trons cooling implies that the electrons temperature 9ei 
is constant along the jet, 9ei{x) ~ 0ei,o- The peak and 
break frequencies thus decay in accordance to Eq. [201 

The observed fiux at frequency v < Vpeakfl can be cal- 
culated by integrating the emitted flux from the differ- 
ent jet segments (this is illustrated in Fig. In order 
to carry an analytical calculation, wc approximate emis- 
sion from a jet segment at position x..x + dx as con- 
tributing only at frequencies v < Vpeakix) (at higher 
frequencies the flux decays exponentially). The decay 
of the magnetic field implies that the peak frequency 
i^peak{x) decays along the jet. Therefore, at a given 
frequency v < Vpeakfij the contribution to the fiux is 
from jet regions only up to a position Xmax for which 
i^peak{xma.x) = V. The dccay law of the magnetic field 
thus gives x^ax = xolv/vpeakfl]"^^""'"' ■ Using the contri- 
bution to the flux from the optically thin regions of the 
spectrum (Eq. [TT|) one obtains 



/3\l/3 2rl uoq^BoXQ ( v \ 

\A) 93^ (l-2ajet/3)meC^mpC^ \Vp„ak,o ) 



dx 



(31) 

The lower integration boundary is taken here as Xmin = 
xo, however the exact boundary is irrelevant as long as 
a^max 3> Xmin- We thus coucludc that for weak magnetic 
field, the observed fiux dX v < Vpeak,o depends on the jet 
geometry via F^, oc (see Fig. [2]). 
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Fig. 2. — Example of the spectra obtained from an initial 
Maxwellian distribution of electrons radiating in weak magnetic 
field, Bq = 3 X 10^ G < Bcr,\, when adiabatic energy losses are 
neglected. The values of the free model parameters are the same 
as the "canonical" values taken in i|2.1.1l Thus, we consider flow 
parameters resulting in no = li ^0 = 45r3, rg = lOrs, and £^,0 = 1- 
We further consider a conical jet, ajet = 1. As the distance to the 
object we take the distance to XTE J1118-I-480, d = 1.8 kpc. The 
inner panel shows the electrons energy distribution at increasing 
distances along the jet. As the electrons temperature is nearly con- 
stant, their density decays as n^j f^t o<: . In the main panel we 
show the resulting flux (solid line), and its decomposition into the 
flux emitted from the different jet segments (dash lines). Marked 
in arrows are the characteristic frequencies Vp^^kfi ^'nd ^thickfi- 
These frequencies are clearly seen when looking at emission from 
an arbitrary jet segments at x > xq, however the overall spec- 
tra below Vpeakfi is flat, Fi, oc (for aj^t = 1, see text). The 
dash-dotted arrows show the direction of the evolution of the elec- 
trons distribution and the emitted flux at different times, which 
are equivalent to distance along the jet. 

For conical jet, a^et = 1, we therefore find that the fiux 
below i'peak,o is constant (i.e., a fiat spectrum, F^ oc j^°). 
W hile this result is siin il ar to the result first obtained 
by iBlandford &: Ko nigl ' (1979") , here it has a different 
physical origin. In Blandford & Konigl (1979) the ori- 
gin of the flat radio spectrum is the decay of the break 
frequency vthick along the jet, while here it is the decay 
of peak frequency, Vpeak ■ The decay law of the break fre- 
quency is similar to that of the peak frequency, and thus 
its existence cannot be observed when the integrated flux 
along the jet is considered. 

We further point that in this regime of Bq < Bcr,i, 
the total (integrated) observed flux below Vpeak,o linearly 
depends on Bq, and on the position of the jet base, xq- 
We show in flgure [2] the results of an exact numerical 
calculation of the observed flux and its decomposition to 
emission from different jet segments obtained in this case. 
The details of the numerical code are given in appendix 

m 

4.2. Intermediate magnetic field: spectral slope 
Fi, cx v'^/"^ at the Optic/UV hands 

In this section we consider intermediate magnetic field 
at the jet base, Bcr,i < Bq < Bcr,2- By definition, for 
Bq > Bcr.i the electrons rapidly cool close to the jet base 
(see §3.11 Eqs. [T5|), while the requirement Bq < Bcr.2 
implies that at the end of the rapid cooling phase, still 
i^thick < ^peak (sBe Eq. I19p . Thc rapid cooling introduces 
an additional break frequency in the observed spectrum, 
i^/ast (see ij3.3I Ea. and figure [3]). 

As explained in and in this case of initial 
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rapid cooling, the analytical calculations are carried in 
two steps. We first calculate the emissivity during the 
rapid cooling phase, under the assumption that the mag- 
netic field is constant during this phase, and is equal to 
Bq. Emission at this phase contributes to the flux at 
the frequency range i^fast < v < Vpeakfi- At a second 
step, we calculate the emissivity during the rest of the 
electrons propagation along the jet, assuming that their 
temperature is constant in this phase (given by Eq. I17p . 
and that only the magnetic field decays. At this stage, 
the electrons emission contribute to the flux at lower fre- 
quencies v < Vfast- The results of the full numerical 
calculations presented in flgure [3l confirms the validity 
of the analytical approximations. 

The fiux at Vpeak.o is estimated by assuming that the 
entire electrons distribution is concentrated at a sin- 
gle energy (characterized by Lorentz factor 7min) and 
that the emission in monochromatic, at Vpeakfi- Thus, 
the total power emitted at Vpeakfi is ^i^|i/=i/peafe o — 
Po/vpeak,Q = 16TTq^BQ/27meC^A. Here, Pq = 
4g^_Bo(7/3)min^/9mgC^ is the total power emitted, and 
the factor A « 10 is introduced here to account for 
the distribution in the electrons energy and for the fact 
that the emission is not monochromatic. The emis- 
sivity is j^\„=,yp,^k,Q = nei,tot X ^V|i/=i/p,„t3.o/47r where 
nei,tot ~ uo/nipC^ is the number density of electrons at 
the base of the jet. The observed fiux from a segment of 
length Ax close to the base of the jet is therefore 



Uo 



(32) 



Comparison with Eq. [3T] derived in §4.11 indeed shows 
that the results are similar up to a numerical factor. 
However, the estimate here is based on robust arguments, 
that can be (and will be) generalized in studying the fiux 
in the power law scenario discussed below. 

The length Ax at which electrons radiation contributes 
to the fiux at Vpeakfi can be estimated by the character- 
istic time it takes the electrons to lose their energy, Eq. 
1141 Since the energy loss is very rapid and occurs within 
time t >tQ, one can approximate the energy loss time in 
the following way. Defining e = t/to — 1 = x/xq — 1, one 
finds that for e <C 1, Eq. [14] can be written in the form 



1 + 



arBltoe 



where the last equality holds for e ^ 
6Trmec/aTBQ^(3{to)to^ The characteristic energy 
loss time of electrons with initial momentum 7/3(to) can 
therefore be written as 

(7tB57/3(*o) 

Writing Ax — jj/3jcAt and using this result in Eq. [32l 
one finds that the fiux at Vpeak.o is 



4(j3(jjAfdisk 



27 (P- nip CUT BoOeiflA. 



(35) 



At frequencies below Vpeakfl and above Vfast, the fiux 
can be calculated as follows. At this frequency range. 

Note that for Bq > -Bci-,ii by definition 
67rmec/oT_BQ7/3(to)to < Ij ^^nd tlius the solution is valid. 



the fiux is dominated by emission from electrons during 
their initial rapid cooling phase. As the electrons cool, 
the peak of the synchrotron emission v = I'peak (t) decays. 
The power emitted at v, -PjyUp^„^(t) is independent on the 
electrons energy (see discussion and Eq. [321 note that 
since B = Bq is constant, Eq. [32] is in fact valid for emis- 
sion at all the frequencies in range Vfast < v < Vpeak.a, 
when the appropriate length Ax is taken). Therefore, 
the observed flux at these frequencies F^ivfast < f < 
Vpeakfi) depends only on the electrons cooling time, or 
cooling length Ax, during which synchrotron emission 
contributes to the fiux at frequency v. Using the de- 
cay law of the electrons energy at early times derived 
above, Eq. [33l one finds 7/3(t) oc {toe)~^ cx Ax~^. The 
peak emission frequency depends on the electrons mo- 
mentum via V — Vpeak{t) oc B{jf3)'^{t) (X Ax~^. We 
therefore conclude that i'{t) oc Ax^^. Using again Eq. 
1321 one immediately concludes that since Fi, oc Ax, then 
for i^fast '£ V < Vpeakfi thc flux decays as Fi, oc 
Note that this result is independent on the jet geometry. 

Once the electrons cool to their asymptotic tempera- 
ture, 9elJ^ they continue to propagate along the jet with- 
out further significant energy losses. For intermediate 
value of the magnetic field, Bcr.i ^ Bq < -Bcr.2, in spite 
of their cooling still vthick < i^peak- As a result, at low 
frequencies v < Vfast, the analysis of the fiux carried in 
§4. II holds. In particular, we conclude that for v < Vfast, 



Fy (X u 



l-l/Ojct 



, and that the fiux at these frequencies 



linearly depends on Bq and xq.^ 

Examples of the spectrum are presented in Fig. [3l 
In the inner inset we show the electrons energy distribu- 
tion in different regions along the jet (the arrows indicate 
the direction of the temporal or equivalently the spatial 
evolution). In spite of the early, rapid cooling, still the 
electrons maintain a quasi-Maxwellian distribution along 
the jet, with asymptotic temperature given by Eq. [TTl 
The main panel shows the observed flux, and its decom- 
position into the flux emitted from the different segments 
along the jet. The flux decay F^, oc at the range 

i^fast < v < Vpeakfi is pronounced. As can be seen by 
the decomposition, contribution to the flux at these fre- 
quencies is only from the innermost parts of the jet. 

4.3. Strong magnetic field: suppression of the fiux at 
radio frequencies 

If the magnetic field at the jet base is higher than i?cr,2, 
then the electrons rapidly cool to temperature below the 
critical temperature 9cr (see Eq. [T8| . Once this happens, 
the peak of the synchrotron emission is in the optically 
thick region, and is therefore obscured, Vpeak < ^tMck- 
The observed break frequency Ujast in this case is given 
by emission from electrons at Ocr and was calculated in 
gS31 Eq. [271 

Emission at lower frequencies that occur as the elec- 
trons continuously cool below 9cr is obscured, since it is 
in regions of high optical depth {vp^ak < vthick)- As a 
result, for Bq > Bcr,2 the peak of the observed fiux is at 
frequency Vfast given by Eq. I27i while at lower frcqucn- 

^ The linear dependence of the flux below vjast on Bq can be 
derived indirectly, from the analysis carried above. It is an immedi- 
ate consequence of the facts that Vpeakfi Bq, Fv\vp^^k g °^ ^0 ' 
Ufast oc Bg''^ and F„\,^j^^^<„^„^^^^^ ^ <x i^-^/^. 
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Fig. 3. — Example of the spectra obtained from an initial 
Maxwellian distribution of electrons radiating in an intermediate 
magnetic field, Bcr.i < Bq = 10^ G ~ Bcr.2, when adiabatic en- 
ergy losses are neglected. Apart from the magnetic field, all the 
other parameters values are the same as in Fig. [2] Inner panel: 
electrons energy distribution at different regions in the jet. The 
dash-dotted arrows show the temporal evolution: at first the elec- 
trons rapidly cool close to the jet base, while maintaining an ap- 
proximate Maxwellian distribution. At a later stage, the electrons 
number density decreases as n^; tot r~^. Main panel: the re- 
sulting spectrum (solid line) and a decomposition to the spectrum 
emitted from different segments (dash lines). We mark the tran- 
sition frequencies Vpeakfl a-nd ^fast- The initial rapid electrons 
cooling results in a fast decay of Vp^aki^), which leads to a fiux 
decay Fv oc at the range Ufast < i' < Vpeakfi- Following 

the rapid decay, the two characteristic frequencies Vthick a-nd ^peak 
have nearly the same values (see text) and evolve in a similar way 
along the jet, as is seen by the decomposition. As a result, below 
'^fasty the flux decays in a similar way to the decay in the weak 
magnetic field scenario, described in Fig. [2] 

cies the flux decays. We therefore conclude that further 
increase of the magnetic field above Bcr,2 results in an 
increase of the observed spectral break frequency, I'fast, 
in accordance to the dependence given in Eq. [27] (as op- 
posed to the result obtained for Bcr,i < Bq < Bcr,2, see 
equation 

Similar to the scenario presented in i i4.2[ in strong 
magnetic field the electrons lose most of their energy 
rapidly. Thus, the calculation of the flux at Vpeak,o in 
Eq. [35] holds. Similarly, below Vpeakfi and above i^fast, 
the flux decays as F„ cx 

At lower frequencies, radio - IR, emission is dominated 
by electrons propagating along the jet with temperature 
0eij < 0cr- Emission from these electrons is charac- 
terized by Vpeak(x) < Vthick{x) . For a Maxwellian dis- 
tribution of electrons, above Vpeak the flux decays ex- 
ponentially. Thus, an exact calculation of vthick re- 
quires s olving a transcend e ntal e quation. Using the re- 
sults of I Jones fc Hardee I ([1979[ ) for z 3> 1 (see dis- 
cussion below Eq. [3]), and the optical depth given 
in Eq. [5] the break frequency vthick can be calcu- 
lated by solving the Eq. [^vthickix) /2vpeak{x)Y^''^ = 
{2/'i)\og['K'^r{x)Aq^c/^vthick{x)]- Here, we approximate 
the logarithm on the right hand side as constant, which 
enables us to write vthick (x) — avp^ak [x) . The parameter 
a depends on Bq . We find numerically that a has charac- 
teristic value of few tens (a « 3 for ~ e_B,cr,2, a w 60 
for tB — 10eB,cr,2 and a w 100 for €b = 100eB^cr,2)- 

With these approximation, one can calculate the ex- 
pected flux at the radio frequencies by integrating the 
flux in the optically thick part of the spectrum (Eq. [T0| 



over the different jet segments, in a similar way to the 
calculation of the radio flux carried in ij4.1l Since the 
electrons cool to their terminal temperature Oeij close 
to the jet base, at a given frequency v < Vthick(x) emis- 
sion is obtained from jet regions xq < x < Xmax, where 
a^max — Xi){i' / avpeakj)~'^^°'^'^ ■ Integrating the flux in the 
optically thick part of the spectrum, using Eq. [TU] one 
obtains 

^ ^m^e^^j:.^ dx 

(36) 

We therefore find that for Bq > Bcr,2, the radio flux 
decreases with the increase of Bq, as Fi, cx O^i jBq (x 

Bq^{(x Cg*). We find numerically that the decay law is 
somewhat weaker, due to the non-linear dependence of 
a on Bq. We further find that the radio flux depends 
on the jet size as Fi, oc Xq"^, and on the jet geometry 
in a similar way as in the weak magnetic field scenario, 
F^ cx jyi-i/%-et. 

Example of the spectra obtained is presented in figure 
[4] Clearly, the flux at high frequencies, above Vfast is 
similar to the flux obtained for weaker magnetic fleld 
(see flgure[3]). However, at lower frequencies, mainly at 
the radio band, the flux is strongly suppressed due to 
the self absorption. This result can naturally explain 
differences in radio emission that are not accompanied 
by similar differences at higher bands in various sources, 
by an adjustment of a single parameter- the magnetic 
fleld. 

In Fig. [5] we compare the resulting spectra for different 
values of the magnetic fleld. While this flgure is derived 
for power law distribution of the energetic electrons (see 
|J5]), the resulting flux below Vpeakfi are independent on 
the distribution of the electrons above 7min and are thus 
similar for both cases. This flgure illustrates how an in- 
crease of the magnetic field from low values first leads 
to an increases in the radio flux, which changes to a de- 
crease in the observed flux when Bq > Bcr,2- The flux at 
the X band on the other hand is much less sensitive to 
the value of the magnetic fleld. We further discuss the 
implications of this in ^ 

5. EMISSION FROM POWER LAW DISTRIBUTION OF 
ELECTRONS WITHOUT ADIABATIC COOLING 

We next consider an initial distribution of the accel- 
erated electrons in the form of a power law energy dis- 
tribution above the peak of the Maxwellian and below 
a maximum Lorentz factor 7max, with power law index 
p. Equation [14] that describes the decay of the electrons 
momentum due to synchrotron emission is independent 
on the electrons distribution, and hence is valid in this 
scenario as well. As a result, most of the analysis carried 
in 2] is valid here. However, the inclusion of a power law 
implies that there is an additional critical value for the 
magnetic field, Bcr.Q (see equation [T6l) . Thus, in cal- 
culating the resulting spectrum one needs to discriminate 
here also between the cases Bq > {<)Bcr,Q, since these 
two cases are qualitatively different. For weak magnetic 
field, Bq < Bcr.Q the electrons distribution in unaffected 
by synchrotron cooling. In contrast, for Bq > Bcrfl, the 
high energy tail of the electrons distribution is cooled sig- 
nificantly, while the low energy tail may be less affected. 
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Fig. 4. — Example of the spectra obtained from an initial 
Maxwellian distribution of electrons radiating in strong magnetic 
field, Bo = 3 X 10^ G > Bcr,2, when adiabatic energy losses are 
neglected. All the other parameters are the same as in Fig. |2] 
Clearly, the flux at high frequencies, v > Vfast is similar to the 
case of lower magnetic field (see Fig. [Sj, while at low frequencies, 
below Vfaat the flux is strongly suppressed. The decomposition of 
the flux to emission from the different jet segments, reveals that 
the origin of this suppression is the suppression of emission from 
the different segments, which are in the optically thick regime fol- 
lowing the initial rapid cooling. The dash-dotted arrow shows the 
temporal evolution of the emitted flux along the jet. Clearly, emis- 
sion at high frequencies is from electrons during their initial fast 
cooling phase, and therefore the only effect of the strong magnetic 
field is a shift in Vfast (see text for details). The electrons temporal 
evolution is very similar to the one presented by the inner panel in 
figure [S] and is thus omitted. 

As a result, the initial electrons energy distribution is 
modified. 

5.1. Weak magnetic field 

We first consider a very weak magnetic field at the 
jet base, Bq < Bcr,o- In this case, electrons cooling in 
the entire energy range is insignificant. T hus, standard 
synchrotron theory (e.g.. lRvbicki fc Lightm an 1979) im- 
plies that above Vpeak and below t'max, the power emit- 
ted by a jet segment is Pi, oc iy~'-P~^^^^. The observed 
flux is calculated in a very similar way to the calcula- 
tion of the flux done in fJH by integrating the contribu- 
tion to the flux from the different segments along the 
jet. In particular, at low frequencies v < i>peak,o, a sim- 
ilar analysis to the one carried in N4.1I holds, and thus 
Fy (X B{){v / Vpeak oY~^^ ""^"^ ■ This is illustrated in figure 

El 

The arguments that led to Eq. [321 can be generalized to 
emission at high frequencies, v > Vpeak from an arbitrary 
jet segment at position x..x -\- dx, as long as electrons 
in this jet segment contribute to the emission at this 
frequency. One therefore finds that the contribution to 
the flux at high frequencies from a jet segment is given 

by 



2q^B gj A/disk 
27T:d'^meC^AjjPjmpC \ vpeak 



-(p-l)/2 



dx. 



(37) 

In order to calculate the total flux, one needs to inte- 
grate over the jet segments that contribute to the flux 
at a given frequency v. Since in this case of weak mag- 
netic field the electrons energy is constant along the jet, 
the peak frequency decay is attributed only to the decay 
of the magnetic field, and Vpeak = i^peakfiix/xoy^^'* . 



The observed fiux at frequency v > Vpeakfl is therefore 
calculated by integrating the emission from jet regions 

(p-i) 



a,et 2 ^ 



1 - 



(38) 

For a given frequency v <C J^max, contribution to the fiux 
is from regions up to Xn^^^{v)/xQ = iy j v^^.^'^ ^ > 1. 
The observed fiux at high frequencies thus depend on 
the jet geometry. For wide jets, a^^t > 2/(]j+ 1), the 
second term in the parenthesis of equation [351 can be 
neglected. The fiux above Vpeak.a is therefore propor- 
tional to Fn cx For narrower jets, 
1/2 < a jet < 2/(p+ 1) and v J^max, the second term in 
the parenthesis of equation 1381 dominates, and one finds 

Fy cx Thecaseoffljet < 1/2 is treated 

in gll 

5.2. Intermediate magnetic field: gradual change of the 
spectral slope at the X band 

If the magnetic field is in the range Ber.o < Bq < Ber,i 
then electrons at 7niax cool significantly at t ^ to, while 
electrons at the lower end of the energy distribution, 7 > 
7min do not significantly cool. This introduces a break 
at vfast, given by Eq. [26l We note that for Bq < Bcr,i, 
Vfast > Vpeakfl- At the range Vpeakfi < V < vfast the 
analysis carried in ^5.11 is valid. Hence, for wide jets 
a jet > 2/{p+ 1), Fy oc Similarly, at lower 

frequencies, v < Vpeak,o, oc v'^^'^/°-i't . 

At higher frequencies, v > v fasti the spectrum is mod- 
ified. Since the energetic electrons rapidly cool, emission 
at these high frequencies occurs only close to the jet base, 
where the magnetic field is constant. The flux at these 
frequencies can therefore be estimated by integrating the 
fiux emitted from a jet segment, given by Eq. 1371 while 
keeping B = Bq and Vpeak = Vpeakfi- 

The length Aa; = x^s,y^{v) — Xq during which elec- 
trons emit at frequency v is related to the cooling rate 
of the energetic electrons. As this cooling rate is in- 
dependent on the electrons initial energy, we can use 
the result derived in §4.21 Eq- Ell to write Aa; = 
jjPjcSt w 67rmeC^7j/3j/crT-So(7/^)max- As the electrons 
cool the maximum emission frequency decays. Since 
the emission frequency v oc (7/?)^, emission at fre- 
quency v can be obtained only from energetic electrons 

(7/3)max(a;)/(7/3)min = / '^peakfi)^^'^ ■ Usiug thcSC rC- 

sults in Eq. [371 one finds 



F, 



ISnd'^mpcA Boq 



'eLO 



-p/2 



(39) 



^ ^peako , 

We therefore conclude that for intermediate magnetic 
field Bcr.o < Bq < Bcr,i, the fiux above Vpeakfi changes 
from Fi, oc i/~(p^i'/2 below Vfast, to oc v~'pI'^ at higher 
frequencies. In reality, we do not expect a sharp cutoff 
at V fasti but a smooth connection (see Fig. [5]). This 
can be the source of the steep spectral slope observed 
in XTE Jl 18+480 at the X-band, and can affect our 
interpretation of the power law index of the accelerated 
electrons. See further discussion in SjQl 
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Fig. 5. — Examples of spectra obtained from an initial power law 
distribution of electrons with power law index p = 2.5, for various 
values of the magnetic field, around the critical values Bcr,i and 
Bcr,2- Bq = 3 X 10^ G (solid, blue), 10* G (dash-dotted, green), 
10^ G (dash-dash, red), and 3 X 10^ G (Dotted, purple). Adiabatic 
energy losses are neglected here. All the other parameters are the 
same as in Fig. [2] The key result here is that an increase of the 
magnetic field above Bcr,2 decreases the radio spectra, with only 
minor effect on the spectra at the X- ray band. Marked are the 
transition frequencies Up^akfi ^^'^ ^Jast- For ^0 > Bcr,2, i^fast in- 
creases with Bq (see Eq. \27h The dashed arrows indicate the evo- 
lution of the spectra with the increase of the magnetic field. At low 
frequencies, below I'peafe.Oi the spectra obtained from Maxwellian 
distribution (Q for similar values of the magnetic field is similar 
to that obtained here. At higher frequencies, above Up^^kfi the 
spectral slope gradually changes from (p — l)/2 for Bq < Bcrfi to 
p/2 for Bo > Bcr.l. 

5.3. Strong magnetic field 

By definition, if Bq > Bcr,i then electrons at the peak 
of the Maxwelhan, hence aU the electrons above the peak 
rapidly cool by synchrotron emission close to the jet base. 
The results derive in the previous sections, M.21 i )4.3l and 
^5.2\ enable us to obtain the emission in the entire spec- 
tral range without the need for additional calculations. 

Since i>fast < i^peakfi, the flux at Vpeakfl is given by 
Eq. [35l while the flux at higher frequencies drop as oc 
At lower frequencies, Vfast < v < Vpeakfi, the flux 
is Fn oc v~^/'^. At even lower frequencies, v < Vfast, the 
flux depends on the jet geometry, F^, oc . For 

intermediate magnetic field, Eq < Bcr.2 the results of 
i j4.2l show that at low frequencies F^, (x Bq, while if the 
magnetic field is strong, Bq > Bcr,2 then the flux at low 
frequencies is suppressed by self absorption, F^, oc Bq^ 
(see MM- 

Comparison of the flux for different values of the mag- 
netic field is presented in Fig. [51 As discussed above, 
regardless of the value of the magnetic field, the low en- 
ergy part of the spectrum, below Vpeak,o is similar in 
this scenario and in Maxwellian distribution scenario dis- 
cussed in SjH We thus conclude that it is not possible 
to discriminate between these two initial distributions 
by observations at low energies, radio - IR bands. At 
higher frequencies, above fpeakfi, the spectral slope grad- 
ually changes from (p — l)/2 for < Bcr,o to p/2 for 
Bq > Bcr.i- However, as is seen in the figure, this tran- 
sition is gradual. We therefore expect an accurate mea- 
surement of the power law index p to be difficult. 

6. MAXWELLIAN DISTRIBUTION OF ELECTRONS WITH 
ADIABATIC ENERGY LOSSES 



While in the previous sections we neglected adiabatic 
energy losses, in this section we consider these, in addi- 
tion to the synchrotron losses. As a result, the electrons 
temperature do not asymptote to a constant value, but 
continuously decreases as the electrons propagate along 
the jet (see Eq. fTS]) . This introduces an additional break, 
t^trans (sec WS.21 equation . 

As we will show below, a major part of the analysis 
carried in the previous sections holds here as well. This 
is due to the fact that in strong magnetic field, the initial 
rapid cooling of the electrons close to the jet base, which 
determines the flux at high energies, is not affected by 
the inclusion of adiabatic losses (see discussion in H2.'3\i . 
These only affect the late time decay of the electrons 
temperature, hence the emission at low (radio) frequen- 
cies. 

6.1. Weak magnetic field: spectral break at the 
transition frequency 

If Bq < Bcr,i, then the initial rapid cooling is insignif- 
icant, and the electrons temperature vary along the jet 
as 9ei oc a;-2ajet/3 ^ggg The magnetic field 

also decays, and as a result, the peak frequency decays 



as I'peak OC X 
Vthick OC X"^" 



7£ijet/3 g^j^j ^j^g break frequency decays as 
"=*/3 (see Eq. Assuming that at 



the jet base Vpeak,o > vtMck.o, the difference in the de- 
cay laws introduces an additional break in the observed 
spectrum, at frequency vtrans (^ ^3.21 Eq. 1^5)) . 

At the frequency range vtrans < v < Vpeakfi, one can 
estimate the observed flux in a very similar way to the 
calculations done in ij4.11 by integrating the emission 
from the different jet segments. The decay of the peak 
frequency implies that contribution to the flux at a given 
frequency v < Vpeak.o is only from jet regions x < Xmax, 
where Xmax/xo = {'^/'^peak,o)~^^''''°'''"^ ■ Integrating the 
emission from the jet segments xq < x < Xmax while 
considering only the optically thin part of the spectrum 
(see Eq. [TT]) gives 



Fy{vtr 



<V<1>, 



'peak^O 



9d2 



(^) 



(I) 



1/3 



1/3 



dx 



9d^ 



(l-2ajet/9)i 
(3/7)(l-l/aje 



(40) 

We therefore find that at the frequency range vtrans < 
< Vpeakfi^ the flux depends on the jet geometry as 
Fi, oc i/(3/'')(i-i/°i=t). This dependence is different than 
in the pure synchrotron case considered in ^14.11 however 
we note that for conical jets {ajet = 1), still a flat spec- 
trum is obtained. At frequencies higher than Vpeakfi the 
flux decays exponentially for a Maxwellian distribution 
of electrons. 

The main difference between this scenario and the pure 
synchrotron scenario discussed in >j4.1l is the appearance 
of the transition frequency, vtrans- This transition fre- 
quency is unavoidable, as it results from the different 
decay laws of Vpeak and vthick that take place as the 
electrons propagate along the jet. By definition, at fre- 
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Fig. 6. — Example of the spectra obtained from an initial 
Maxwellian distribution of electrons radiating in weak magnetic 
field, Bq = 3 X 10^ G < Bcr,i, when adiabatic energy losses are 
considered. All the other parameters are the same as in figure [2] 
The dash-dotted arrows show the direction of the evolution of the 
electrons distribution (inner inset) and flux (main panel). The elec- 
trons main cooling mechanism is adiabatic energy losses, which re- 
sults in 9^1 oc The transition frequency utrans is clearly 
marked. The decomposition into emission from jet segments, il- 
lustrates that the transition occurs once the emission from th e je t 
segments becomes optically thick, above xtrans (see equation I22I I. 
At Utrans <U < Up^ak.O, F„ OC u^^/'''>(^-^ / "-j't) ^ which for ajet = 1 
as considered here results in flat spectra. However, below utrans a 
flat spectra cannot be obtained (see text for details). 

quencies v < vtrans, the peak of the emission is in the 
optically thick part of the spectrum {vpeak < ^thick), and 
is thus obscured. The flux at low frequencies v < vtrmis 
is therefore dominated by emission at vthick- This sce- 
nario is similar to the one discussed in iJ4.3l We showed 
there that for v^eak < vthick-, the evolution of vthick is de- 
termined by solving the Eq. [^i'thick{x)/2vpeak{x)Y^^ = 
(2/3) log[7r^r(a;)yl(7^c/4z/tfticfc(a:)]. A rough estimate of 
the flux at low energies can be done by approximat- 
ing the logarithm as constant, which enables to write 
vthick ot Vpeak OC x~^°J=*/'^. Sincc in this case of emis- 
sion in the optically thick part of the spectrum, dFn oc 
v^rOeidx cx v'^r^^^dx (see Eq. [TU)) . we conclude that 

An example of the obtained spectrum is illustrated in 
figure [6] The decomposition shows the evolution of the 
characteristic frequencies vthick and Vpeak along the jet. 
The faster decay of Vpeak is clearly pronounced. We mark 
the transition frequency vtrans- The different spectral 
regimes are easily identified in this plot. 

6.2. Intermediate and strong magnetic field 

For higher values of the magnetic field > Bcr.i, the 
electrons rapidly cool by synchrotron emission close to 
the jet base. Only after the initial rapid cooling, the 
electrons temperature decay law asymptotes to 9ei oc 
^-2aj,t/3^ During the initial rapid decay the adiabatic 
energy losses are negligible (see H2.31 equation [13]) . This 
enables us to use the analysis carried in H4.2\ and M.3\ in 
calculating the flux at high frequencies, which are dom- 
inated by emission during the initial rapid decay phase. 
On the other hand, the flux at low frequencies is governed 
by emission from electrons as they propagate along the 
jet. This enables us to use the analysis carried for the 
case of weak magnetic field in H6.1\ in calculating the flux 
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Fig. 7. — Example of the spectra obtained from an initial 
Maxwellian distribution of electrons radiating in intermediate mag- 
netic fleld, Bcr.l ^ Bo = 3 X 10* G < Bcr,2, when adiabatic energy 
losses are considered. All the other parameters are the same as 
in Fig. [2] The dash-dotted arrows show the direction of the evo- 
lution of the electrons distribution (inner inset) and flux (main 
panel). The initial rapid decay of the electrons temperature is seen 
in the inner panel. While we mark in this figure the three transi- 
tion frequencies, u^^^i^ q, Uf^^^t and utrans, clearly for intermediate 
value of Bo the transition in the flux behaviour is smooth. 

at low frequencies. 

For intermediate magnetic field Bcr.i < Bq < Bcr.2, 
similar analysis to the one carried in i)4.2l shows the ex- 
istence of a break frequency at Vfast, whos value is given 
in Eq. [261 At the frequency range vfast < v < Vpeak.,o, 
the flux decays as Fy oc and at higher frequencies, 

V > Vpeakfi, the flux decays exponentially. 

The initial rapid cooling results in a shift of vtrans 
to higher frequency, given by Eq. [24l As discussed in 

i i3.3l for Bq < Bcr,2, Vtrans < Vfast- SiuCC bcloW Vfast 

the analysis carried in holds, we conclude that at 

the range Vtrans < V < Vfast, cx i/-(V3)(l-l/ajet)^ 

while at lower frequencies f < vtrans, F,j{v < vtrans) oc 

j,(13-3/a,-et)/7^ 

We thus conclude that for Bcr.i < Bq < Bcr,2, the 
spectra below Vpeak,o has three different regimes, sepa- 
rated by the two transition frequencies vtrans and Vfast- 
In practice, however, we expect these transitions to be 
smooth, as is illustrated in Fig. [7] 

If the magnetic fleld is stronger, Bq > Bcr.2, then, by 
deflnition emission below i^fast is in the optically thick 
part of the spectrum. In this case, vtrans = Vfast (see 
discussion in H3.3p . Thus, the spectrum below fpeafc.o in 
this case is composed of only two segments: F^, cx j/^^/^ 
for Vfast < V < Vpeak,o, and Fy cx iy(i3-3/a3et)/7 lower 
frequencies. This is illustrated in Fig. [5| We note that 
in this scenario, the requirement for a flat radio spectra 
(F^ cx below Vtrans) Can not be fulfilled for aj^t > 1/2. 

7. POWER LAW DISTRIBUTION OF ELECTRONS WITH 
ADIABATIC ENERGY LOSSES 

When adiabatic energy losses are included, the elec- 
trons power law distribution above the Maxwellian af- 
fects not only the high energy part of the spectrum {v > 
Vpeak,a), but also the very low energy part, v < vtrans as 
well. This is due to the inclusion of vio^, as discussed in 
El 

The analysis of the spectra in this case is very similar 
to the analysis already done in the former sections. The 
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Fig. 8. — Example of spectras obtained from an initial 
Maxwellian distribution of electrons when adiabatic energy losses 
are included, for various values of the magnetic field, around the 
critical values Bcr,i and Bcr,2- Bq = 3 x 10^ G (solid, blue), 
3 X 10* G (dash-dotted, orange), 10^ G (dash-dash, red), and 
3 X 10^ G (dotted, purple). All the other parameters are the same 
as in Fig. [2] As the magnetic field increases, vtrans shifts to higher 
frequencies, and for Bo > Bcr,2, vtrans = vjast- As a result, a fiat 
radio spectra cannot be obtained at low frequencies in this model. 
The increase of the magnetic field resu lts in a shift of Vfast to 
higher frequencies, as predicted by Eq. 1271 For Bq > Bcr,i, the 
flux at higher frequencies is nearly insensitive to the exact value of 
the magnetic field. 

observed spectra is composed of several break frequen- 
cies, in accordance to the value of the magnetic field. 

For weak magnetic field, Bq < -Bcr,o, the spectra 
is composed of four distinctive segments, separated by 
three break frequencies: fpeafc, Oj I'trans and vi^w, given 
by Eqs. S] , [53] and [251 Above Vpeakfi and below i^max,0: 
the flux oc i/-(p-i)/2 is given by Eq. [SSjin [JSlIl At the 
frequency range vtrans < v < Vpeak.o, the results derived 
in SjHUhold, and oc ,y(3/7)(i-i/a;,0. 

At lower frequencies, viow < v < vtrans, the flux can 
be calculated in the following way. We showed in ij3.3l 
that the spectrum in this range is determined by emis- 
sion from a power law distribution of electrons at the 
break frequency vtMck which marks the transition from 
the optically thin to the optically thick part of the spec- 
trum. In the optically thick part of the spectrum (Eq. 
[TU)) . dF^ oc r'^{j^/Tv)dx oc rB~^^'^v^/'^dx, where jjj and 
are the emissivity and optical depth for power law 
distribution of the emitting electrons (see §3.3|) . Using 
V = vthick, and the dependence of vthick on the position 
along the jet as derived in i^thick oc r^^^ , one finds 

di^^ oc (41) 

where Ai{p) = (8 -f 7p)/{3p+ 12) was derived in ij3.3l 
T he result in E q. [4T] is identical to the result derived 
by iKaiserl ()2006il : however, we note here that this result 
is valid only at the frequency range z^joji, < v < vtrans ■ 
Writing F^, oc at this range, we find that for p ~ 2 
(2.5), for ajet = 1, = 0.45(0.59), for Ojet = 3/4, 
T] = 0.18 (0.33), and for a^et = 2/3, 77 = 0.05 (0.20). One 
can thus conclude that for a jet ^ 2/3, a flat spectrum 
is expected at t his range viow < v < vtrans, similar to 
the conclusion in IKaiserl (|2006l ). At much lower frequen- 
cies V < flow the discussion in §6.21 holds, and thus we 
can approximate oc i/(i3-3/ajet)/7^ high frequen- 
cies, V > Vmaxfi, an exponential decay in the spectrum 
is expected. 



From the discussion in [JSJ we find that the spectrum 
obtained for intermediate value of the magnetic field, 
Bcrfi < Bq < Bcr,i differs than the spectrum obtained 
for Bq < Bcr.i, by the inclusion of Vfast (equation [26]) . 
Since Vfast > t^peak,o, this inclusion only affects the ener- 
getic part of the spectrum: for v > Vfast, oc v~'^p/'^\ 
while the spectrum at lower frequencies is similar to the 
spectrum obtained in the weak magnetic field case de- 
rived above. For magnetic field in this range, the values 
of ^trans and viow are determined by Eqs. 1241 1301 respec- 
tively. 

For higher values of the magnetic field, Bcr,i < Bq < 

Bcr,2, Vfast < l^peakM and l^iow = Vtrans- Thc SpCCtrum 

is thus separated to four distinctive regimes, which are 
different than those for the weak magnetic field derived 
above. For ly > Vpeak,o, Fv ^ v~'^p/'^\ and at lower fre- 
quencies Vfast < V < Vpeak,s, Fi, OC i/^^^^. At cveu lower 
frequencies, vtrans < v < Vfast, F^ oc i/(3/7)(i-i/ij.t)^ 
while at the low end of the spectra, v < vtrans, F^ oc 

j,(13-3/Q,et)/7^ 

For even higher values of the magnetic field, Bq > 
Bcr,2: Vtrans = Vfast are givcu by Eq. [21] For this 
value of the magnetic field, there are only three distinc- 
tive spectral regimes below Vmax.o- v < Vfast, Vfast < 
V < Vpeak,o and v > Vpeak,o- The flux in these regimes is 
similar to the flux in the equivalent regimes obtained for 
weaker magnetic field. 

Examples of the spectra obtained in this scenario are 
presented in Figg. [9land[T0l In Fig. [9]we present the de- 
pendence of the spectra on the strength of the magnetic 
field. We mark thc transition frequencies. Note that for 
strong magnetic field there is a degeneracy: e.g., as ex- 
plained above, for Bq > Bcr,2, Viow = Vtrans = Vfast- In 

figure[TO]we give examples of the spectra obtained for the 
several possible jet geometries. Clearly, for low values of 
a jet, the flux at low frequencies, at the range viow < v < 
Vtrans iucrcascs (this holds only as long as ajet > 1/2), 
while the flux at high frequencies, v > Vpeak.o is unaf- 
fected. The theoretical approximation in Eq. [^predicts 
that for power law index p = 2.5 as is used in the plots, 
a flat spectra is obtained for ajet = 0.56. The numerical 
result is in very good agreement with this prediction. 

8. NARROW JETS 

As discussed in ij2.31 for narrow jets ajet < 1/2, the 
electrons lose their momentum asymptotically as 7/? oc 
j^2ajEt-i^ This asymptotic decay law is independent on 
thc inclusion of adiabatic energy losses. Therefore, the 
analysis is similar in both scenarios. Moreover, due to 
the continuous electrons energy loss along the jet, the 
analysis is essentially similar to the analysis carried in 
fJS] f|7]when adiabatic energy losses were included. We 
give here a brief description of the spectrum obtained in 
this case. As we will show below, we find it impossi- 
ble to obtain a flat radio spectra, as suggested by many 
observations. We thus find this scenario less likely. 

The high energy part of the spectra, at > 
min(t'peafc,o, Vfast) is determined by electrons close to the 
jet base, and is therefore unaffected by the jet geometry. 
The analysis in the previous sections are thus valid in 
this case as well. 

The difference in the flux emitted in this scenario 
and the scenarios discussed above is thus expected at 
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Fig. 9. — Examples of spectra obtained from an initial power law 
distribution of electrons with power law index p = 2.5 when adia- 
batic energy losses are included, for various values of the magnetic 
field, around the critical values Bcr,l and Bcr,2- Bq = 3 x 10"^ G 
(solid, blue), 10* G (dash-dotted, green), 3 X 10* G (dash-dash, 
orange), and 10^ G (dotted, red). All the other parameters are 
the same as in Fig. [2] We mark the transition frequencies I'peafe.Oi 
^faat, I'trans and Vion,- At the frequency range uio^ <v < vtrans, 
the flux dependence on the jet geometry is described in Eq. 1411 At 
other frequencies the spectral shape is similar to the spectral shape 
discussed in former sections. In particular, the steepening of the 
slope at high frequencies from (p — l)/2 to p/2 with the increase 
of -Bo is pronounced. 
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Fig. 10. — Examples of spectra obtained from an initial power 
law distribution of electrons with power law index p = 2.5 when 
adiabatic energy losses are included, for different jet geometries. 
A conical jet aj^t = 1-0 (solid, blue), aj^t =2/3 (dash-dotted, 
green), aj^t = 0.56 (dash-dash, orange) and ttjet = 0.5 (dotted, 
red). The magnetic field is taken as Bq = 3 x 10^ G < Bcr,i, and 
all the other parameters are the same as in Fig. |2] For a^et = 0.56, 
the analytical result in Eq. 1411 predicts a fiat spectrum Fi, oc 
at the range Uiom < ^ < ^trana- The numerical result shows very 
good agreement with this prediction. 

lower frequencies, v < min{iyfastT i^peak,o)- The flux 
at vtrans < V Can be calculated from the decay law 
of the electrons energy. At this frequency range, the 
peak of the synchrotron emission from any jet segment 
is in the optically thin part of the spectrum. As a re- 
sult, one can write z^peafe oc {'-f(3)^B oc x^"-^"*"^, and 
using Eq. [TTl oc r^n{r)Bdx oc x^""^''*. One there- 
fore concludes that at this frequency range Fy oc z/''% 
with 7^1 = (1 — ajet)/{'iajet — 2). For ajet = 0, we find 
771 = —1/2 while for a^et = 1/2, 771 = —1. 

If the electrons have a Maxwellian distribution, then 
at lower frequencies, v < vtrans the peak of the emit- 
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Fig. 11. — Examples of spectra from power law distribution of 
electrons with index p = 2.5 for different jet geometries in narrow 
jets. Shown are ajet = 0.1 (solid, blue), aj^t = 0.2 (dash-dotted, 
green), aj^t = 1/3 (dash-dash, red). The magnetic field is taken as 
Bo = 3x10^ G < Bcr, 1 1 and all the other parameters are the same 
as in Fig. [2] The transition frequency vtrans is clearly seen. In no 
case a fiat radio spectra can be obtained (see text for details). 

ted spectrum is obscured. Using Eq. [TOl one finds 
dFy oc v^rOeidx oc x'^"-^"^'^. We have used here sim- 
ilar arguments to the ones used in ^6.1[ that enables 
us to estimate the x-dependence of vtMck as similar to 
that of Vpeak at low frequencies. We therefore con- 
clude that at this frequency range, Fi, oc u^^ , with 
772 = (9ajet - 4)/(3ajet - 2). For ajet = 0, 7/2 = 2 
while for Ujet = 1/3, 772 = 1. For somewhat wider jet, 
a jet = 4/9, 772 — 0, and a flat spectrum is obtained. 
However, as pointed out in for aj^t ~ 0.45, the 

transition frequency is very low, Vtrans ~ 10^ Hz. As a 
result, the flat part of the spectra occurs at frequencies 
well below current observation capabilities. 

If the electrons have a power law distribution above 
7min, then the observed spectrum below vtrans shows 
an additional break, at viow Similar to the discussion 
in fJTl for V < viow, oc v^'^ . The flux in the in- 
termediate range, vio^ < v < vtrans is calculated as 
follows. An analysis similar to the one carried in H3.3I 
shows that in this case, when writing the electrons dis- 
tribution as nei{'j)d'y = k{x)"f~Pd'y, the proportionality 
constant evolves as k{x) oc resulting 
in i^thick oc a;K=t(3p-8)+2(i-p)]/(p+4)^ For p = 2 this re- 
sults in vthick oc a;-(i+"j<!t)/3^ while for p = 2.5 one ob- 
tains vthick oc a::-(6+aj.t)/i3^ Using dF^ oc rB^^/'^v^/'^dx, 
one finds that for p = 2, dF^ oc jy-(i+4aj.t)/(2+2a,.t)^ 
while for p = 2.5, dF^ oc j,(4-34a,.t)/(i2+2a,,t) . We 
thus conclude, that for p = 2, the flux varies from 
Fi, oc for ajet = 0, to F^, oc for ajet = 1/2. 

For p = 2.5 a flat spectrum can actually be obtained at 
this range viow < v < vtrans, but only for very narrow 
jets, ajet = 0.1. For such a narrow jet, the analysis car- 
ried in §3.21 §3.31 shows that viow S ^trans, and thus the 
flat spectrum is limited to a very narrow band, and in 
practice does not exist. 

Examples of the spectra for the various jet geometries 
are presented in Fig. [11] The transition frequency vtrans 
is clearly seen to evolve to lower frequencies when ajet 
increases, as is predicted in i|3.2l 

9. SUMMARY AND DISCUSSION 
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In this paper, we have extensively studied a model for 
synchrotron emission from jets in black-hole X-ray bina- 
ries. Our basic assumption is that the electrons are accel- 
erated once at the base of the jet, and lose their energy 
by synchrotron emission and possible adiabatic energy 
losses, as they propagate along the jet. As the details 
of the acceleration process are not understood from first 
principles, we considered two scenarios which have strong 
theoretical and observational motivations: a Maxwellian 
distribution of energetic electrons, and a power law dis- 
tribution above the low energy Maxwellian, at the en- 
ergy range 7min < 7 < 7max, with an exponential decay 
at higher energies. The inclusion of a low energy cutoff 
to the accelerated electrons energy distribution directly 
implies a characteristic break frequency in the observed 
spectra, at v = Vpeak,o (see Eq. S]). This break fre- 
quency adds to the break frequency vthick which marks 
the transition from the optically thin to the optically 
thick emission. 

We assume that the magnetic field decays along the jet 
in accordance to Poynting flux conservation, B{r) oc r~^. 
We showed in ^2.3[ that the electrons cooling along the 
jet has an analytical description (see Eqs. [T3l and [T4|) . 
These equations hold the key to the rest of the analysis. 
By studying the cooling rates, we found that in strong 
magnetic field, the cooling can be separated into two dis- 
tinctive regimes: an initial rapid cooling that takes place 
close to the jet base, during which the electrons emit at 
high frequencies (UV, X- and 7-rays), and a secondary 
long phase in which the electrons cooling rate asymp- 
totes. For wide jets, ajet > 1/2, when only synchrotron 
emission is considered, the electrons energy becomes time 
(and x) independent, 7 cx when adiabatic energy 
losses are considered, then 7 oc x"^"^'''/^. For narrow 
jets, both scenarios result in a similar asymptotic decay 
law, 7 oc x'^'^^^*~^. 

These results enabled us to define in fJXl] critical values 
of the magnetic field, and in ^3.2\ , H3.3[ five transition 
frequencies in the observed spectra. We then analyzed 
the resulting spectra for the various assumptions on the 
accelerated electrons spectra. Maxwellian distribution 
was considered in []4] and [JU and power law spectra was 
considered in Sj5]and SjTl We further studied the effect of 
inclusion of adiabatic energy losses in !jS]and 33 and the 
unique scenario of narrow jets, a jet < 1/2, in [J8l 

While the variety of spectra that can be obtained in 
such a simplistic model is found to be very large, we can 
point to some general properties of the spectra, that we 
find of high importance. 

(I) A flat radio spectra, as it is seen in many objects, 
can be obtained if only synchrotron emission is consid- 
ered, for conical jets (see 52J !j5l Figg. [2]-[5]), regardless 
of the value of the magneti c field. This scenario issimi- 
lar to the orig inal model of iBlandford fc Konigl I (fl979l) . 
although there is a different physical origin for the flat 
spectrum: here the flat spectrum results from a decay of 
the peak of th e emission frequency Vpenk along the jet, 
while in Blandford fc Konigl I (119791 ) model, the evolu- 
tion of Vthick was considered. As we showed in H3.21 ne- 
glecting adiabatic energy losses, both frequencies evolve 
in a similar way. We further showed that when adia- 
batic energy losses are included, a flat spectrum can be 
obtained for conical jets only above vtrans, whose value 



is given in Eq. [23] (see fjH g3 Figg. [6]- [9]). At lower 
frequencies, a flat spectrum can be obtained only when 
a series of conditions are met: the electrons are power 
law distributed, the jet has a specific geometry (see ^\ 
equation Hi] and Fig. [TO]) . and the magnetic field is lim- 
ited. In this case, a flat spectrum is obtained only at the 
range viow < v < vtrans- For narrow jets, we showed in 
S|8]that a flat radio spectra cannot be achieved (see Fig. 
E]). 

(II) We showed that the flux at the radio wavelengths 
depends on the value of the magnetic field in a non-trivial 
way: for magnetic field at the jet base Bq < Bcr,2, the 
flux increases with the increase of the magnetic field (see 
tj4.ip . However, for stronger magnetic field, Bq > Bcr,2 
a further increase of the magnetic field at the jet base 
leads to a rapid decrease of the observed radio flux, due 
to an increase in the optical depth (see HA.'d^ . This 
is not accompanied by a similar change in the flux at 
higher frequencies, v > Vfast'- the flux at high frequen- 
cies asymptotes to a constant value in a strong magnetic 
field, B > Bcr,i ■ We therefore find a natural mecha- 
nism that can lead to a variation in the ratio of the radio 
to X ray fluxes, by a simple change in the value of the 
magnetic field. We further investigate the consequences 
of this idea in a forthcoming paper [Casella & Pe'er, in 
preparation]. These results imply that the strongest ra- 
dio emission occurs for Bq w Bcr,2, or about three orders 
of magnitude below equipartition value. While the origin 
of the magnetic field in jets is not understood from first 
principles, we believe that this value could be used as 
a guideline for models of magnetic field production near 
the jet base. 

(III) We showed that for intermediate values of the 
magnetic field Bq > Bcr.o, the flux at X-ray wavelength 
gradually changes from (p— 1)/2 at low energies to p/2 at 
higher energies, due to the rapid cooling of the most ener- 
getic electrons (see tj5.2p . For an even higher value of the 
magnetic field, Bq > Bcr,i, we found that at the range 
t^fast < V < Vpeakfi, which is typically in the optical-to- 
X-ray band, the flux decays as oc i>~^/'^ (see ij4.2p . 
This decay law occurs close to the jet base, and is there- 
fore independent on the jet geometry. 

This gradual transition of the flux implies that every 
attempt to determine the power law index p of the accel- 
erated electrons by fltting X-ray data should be done 
with great care. The result Fi, oc above Vfast 

is very robust, as it is independent on the initial dis- 
tribution of the electrons, or on the exact value of the 
magnetic field. It can easily be misinterpreted as due 
to synchrotron emission from a power law distribution 
of electrons, whose cooling is insignificant. Under this 
interpretation, one would come to the wrong conclusion 
that the power law index p of the accelerated electrons 
is p = 2. Similarly, at higher frequencies, the gradual 
change of the spectral slope from (p — l)/2 to p/2 as 
predicted here, if fitted using a single power law over a 
limited band, can lead to a wrong conclusion about the 
value of the power law index p. This may be the source 
of the discrepancy between measurements of the index p 
in BHXRBs and in other objects, as discussed in [JT] 

(IV) Wo obtained high radio flux, ~ 10 mJy for param- 
eters that characterize emission from XTE J1118-I-480, in 
models of wide jets. This is similar flux to that observed 
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in this object (|Hvnes et al. Il2000l) . While we did not aim 
to fit data in this manuscript, we find these results en- 
couraging. We further point out that the terms "wide" 
and "narrow" jets used here can be somewhat mislead- 
ing, since they refer only to the confinement of the jet. 
A "wide" jet as defined here can be geometrically very 
narrow. 

Our model is of course far from being able to describe 
the full physical processes that are expected to occur 
inside the jets. Several radiative processes, like Comp- 
ton scattering, pair production or the full effect of syn- 
chrotron self absorption on the electrons energy distri- 
bution, are not considered here (see iR)) . However, these 
phenomena have only minor effects on the resulti ng spec- 
tra u nder the conditions assumed here (see, e.g.. iKaiser I 
I2006D . Moreover, we did not consider in this work the 
possible contribution of internal energy dissipations (e.g., 
internal shocks) that can lead to multiple accelerations 
episodes of electrons as they propagate along the jet. 

Multiple acceleration episodes of electrons along the 
jet results in a complicated spectra, whose details de- 
pend on the details of the acceleration processes (e.g., 
position, fraction of particles that are being accelerated, 
strength of the magnetic field at the acceleration sites, 
etc.). The spectra obtained in this work can therefore 
be viewed as a basic ingredient of the spectra that re- 
sults from such multiple accelerations: spectrum in the 
more complicated case can be obtained by a composition 
of the spectra presented here. In this manuscript, we 
focused on the spectra that results from a single acceler- 
ation episode, in order to demonstrate the key physical 



processes that occur in the plasma. While our numerical 
model can very easily be generalized to include multi- 
ple acceleration episodes, due to the expected complex- 
ity of the spectra in this case, we leave this for a future 
work. We do stress though, that any model that consid- 
ers particle acceleration to high energies should treat this 
phenomenon separately than synchrotron and adiabatic 
cooling, since these phenomena have different physical 
origin. 

One of the key uncertainties in models of emission from 
jets, when internal dissipation of energy takes place, in- 
volves the origin of the magnetic field. Here, we assumed 
that the magnetic field originates at the core, and evolves 
according to Poynting flux conservation law. However, 
an alternative scenario may be that the magnetic field 
is produced by the internal shock waves. In such a sce- 
nario, the magnetic field would evolve in a different way, 
thus some of the results derived here would no t hold. 
Such a scenario was considered bv iKaiserl (|2006f ). where 
however a low energy cutoff in the accelerated electrons 
distribution was not included, which has a strong influ- 
ence on the obtained results. A full treatment of the 
spectral dependence on the origin of the magnetic fleld 
is left for future work. 

We discussed in this work emission from jets in the 
low/hard state of BHXRBs candidates, since in this state 
there are good evidence for the existence of jets. How- 
ever, the physics of emission from jets in other sources, 
like AGN's is most probably very similar. We thus expect 
that many of the results found here should be relevant 
for jets in other astronomical sources as well. 



APPENDIX 
NUMERICAL CODE 

The numerical code is based on the code constructed by iPe'er fc WaxmaiTI ()2005f ) in the study of emission from 
GRB's. The particles are assumed to be accelerated at the base of the jet. A fraction < ep/ < 1 is assumed to 
have a power law distribution between 7niin and 7max, while 1 — Cpi of the particles have a low energy Maxwellian 
tail below 7niin- The particles energy distribution is discretized in momentum space (7/?) in order to obtain accurate 
calculations when the electrons cooling is significant. The code gets as an input the parameters of the flow (see ^2.1^ . 
and calculates the number and energy densities of the electrons at the jet base. 

The code divides the jet into segments, and evolves the magnetic field at each segment, in accordance to B{x) oc 
x-°-i^t , Both Bq and a jet are considered as free parameters. Given the strength of the magnetic field, the code 
calculates the rate of electrons cooling as they propagate through each segment (Eq. [T^ . In parallel, the emitted flux 
from each segment (Eq. [SJ is calculated self consistently by solving the radiative transfer eq uation. When calculating 
the emissivity, the full cyclo-synchrotron spectra from each momentum bin is calculated (|Mahadevan et al~l 119961 : 
IPe'er fc Waxma^ [20051) . Since the magnetic field varies along the jet, the emissivity is calculated for a reference 
value of the magnetic field and is tabulated. As the magnetic fleld evolves along the jet, the characteristic frequencies 
decay; however, using the scaling law for the emitted frequencies, v = vq x (B/Bq) and emitted power, P{v) = 
PliyQ X {B / Bo)]{B / Bq) the full emissivity at any given jet segment is readily obtained by interpolating from the saved 
values. Given the e missivity at each segment, the self absorption coefficient is calculated by solving the integro- 
differential equation (|Rvbicki fc Lightman I [19791 eq. 6.50), which is correct for a general distribution of electrons. 
Calculation of the optical depth followes the assumption that the line of sight is perpendicular to the jet axis. 

When adiabatic energy losses are considered, Eq. ^ which is accurate for any value of the electrons momentum is 
used. The discrete electrons momentum levels are kept constant. Therefore, calculation of the electrons cooling is 
done in two steps. The first step is to obtain their new momentum, by solving equation [2] At a second step, the new 
electrons distribution is fltted to the original, discrete cells. We use a second order fit in logarithmic space of both the 
electrons number and their momentum, which is proven to conserve both number and energy density. 

In this version of the code we do not assume any further interactions between the photons and the electrons (e.g., 
Compton scattering, increase of the electrons energy due to synchrotron self absorption, pair production, etc.). We 
therefore found that a simple first order integration scheme is sufficient in calculating the evolution of the electrons 
energy distribution. 
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